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ABSTRACT 


A  detailed  literature  survey  on  two-phase  flow  of  gas-liquid 
mixtures  in  horizontal  pipes  is  presented. 

Flow  characteristics  of  air-water  mixtures  were  studied  in  a 
1.026  inch  inside  diameter  transparent  plastic  pipe.  The  input  air-water 
volume  ratios  were  varied  from  0.1  to  200  in  some  and  up  to  'JOO  in  other 
cases,  for  10  superficial  velocities  of  the  reference  phase  (water)  ranging 
from  0.01  ft.  per  sec.  to  5*03  ft.  per  sec. 

A  correlation  of  the  experimental  pressure  drop  for  the  flow 
rates  and  ratios  investigated  is  presented  in  terms  of  a  superficial 
friction  factor  based  on  the  properties  of  the  reference  phase  and 
reference  phase  Reynolds  number.  The  data  showed  good  agreement  with 
the  Lockhart  and  Martinelli  correlation,  though  in  the  viscous-turbulent 
and  viscous-viscous  mechanisms  of  flow,  the  present  data  fell  below  the 
Lockhart  and  Martinelli  line. 

Some  photographic  and  mostly  visual  observations  were  made  of 
the  flow  patterns  that  existed  during  the  experiments.  A  tentative  chart 
is  presented  that  enables  the  prediction  of  the  flow  pattern  that  would 
exist  in  air-water  systems  at  given  mass  velocities  of  the  two  phases. 

The  hold-up  ratio  is  presented  as  a  function  of  the  superficial 
velocity  of  the  second  phase  at  different  superficial  velocities  of  the 
reference  phase. 

Slip  velocity  is  presented  as  a  function  of  the  input  air-water 
volume  ratio.  At  zero  input  volume  ratio,  that  is  all  liquid  flow,  the 
hold-up  ratio  approaches  unity,  while  the  slip  velocity  approaches  zero. 
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INTRODUCTION 


In  the  last  decade,  there  has  been  a  decided  trend  towards 
the  transport  of  liquids  and  gases  simultaneously  through  pipe  line 
systems,  more  particularly  in  the  case  of  oil  and  natural  gas.  In 
areas  like  swamps,  offshore  locations  and  hilly  regions  where  trans¬ 
portation  is  generally  a  problem,  this  method  of  transportation  is 
more  favourable. 

Two -phase  flow  is  not  confined  to  the  pipe  line  systems  alone 
but  is  encountered  in  many  engineering  processes.  Horizontal  tube 
evaporators,  chemical  reactors,  heat  exchangers,  dry  return  of  a  vapour 
heating  system,  oil-gas  wells,  boiling  water  reactors  and  the  injection 
of  a  liquid  drying  agent  in  natural  gas  lines  to  reduce  hydrate  formation 
are  some  examples  that  involve  two -phase  flow.  In  some  cases  as  in 
chemical  reactors  it  is  often  a  question  of  necessity  rather  than  choice. 

The  design  of  a  pipe  line  or  any  other  equipment  for  two- 
phase  flow  depends  in  large  measure  upon  an  ability  to  predict  the 
energy  losses  that  can  be  expected.  This  is  possible  only  if  the 
mechanism  of  two-phase  flow  is  well  understood.  The  energy  losses  for 
the  flow  of  single-phase  fluids  can  easily  be  determined  by  the  use  of 
the  well  known  friction  factor  -  Reynolds  number  -  relative  roughness 
relationships  and  the  Fanning  equation.  In  the  case  of  two-phase  flow, 
the  presence  of  the  second  phase  increases  the  number  of  variables  and 
changes  the  geometry  of  the  system.  Consequently  the  problem  of  analysis 
is  more  complicated.  As  yet,  no  complete  and  satisfactory  correlation 
is  available  for  the  prediction  of  the  pressure  drop  resulting  under 
a  given  set  of  conditions.  Moreover  the  relationship  of  the  flow  types 
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or  patterns  to  the  relative  flow  rates  of  the  two  phases  and  their 
physical  properties  is  not  well  understood.  This  much  is  certain 
however  -  the  energy  losses  encountered  in  two-phase  flow  are  ordinarily 
greater  than  those  for  single-phase  flow  (3,  29,  34,  4l).  The  increase 
in  energy  losses  is  generally  attribulted  to  two  factors.  One,  the 
irreversible  work  done  by  the  gas  on  the  liquid  as  it  is  dragged  along  the 
pipe,  and  two,  the  smaller  areas  of  flow  available  to  each  phase  as  a 
result  of  the  presence  of  the  other  phase.  Since  in  turbulent  flow  the 
pressure  drop  is  known  to  be  inversely  proportional  to  about  the  fifth 
power  of  the  effective  diameter,  even  a  small  decrease  in  the  flow  areas 
causes  a  substantial  increase  in  the  pressure  drop. 

Previous  Investigators: 

One  of  the  earliest  experimental  investigations  on  two-phase 
flow  was  reported  by  Boelter  and  Kepner  (8)  in  1939*  They  measured 
the  pressure  drop  for  the  horizontal  flow  of  air-water  and  air-oil 
mixtures  in  a  study  of  a  fuel  distribution  system  for  heaters  in  California 
citrus  orchards.  The  test  sections  consisted  of  50  feet  long  1/2  inch 
and  3/4  inch  steel  pipes.  They  suggested  that  the  pressure  gradients 
were  functions  of  several  dimensionless  groups  but  no  correlation 
was  presented. 

Martinelli,  Boelter,  Taylor,  Thomsen  and  Morrin  (29)  presented 
the  first  more  or  less  general  correlation  in  1944,  for  gas  turbulent- 
liquid  turbulent  and  gas  turbulent-liquid  viscous  flow  mechanism.  A 
further  analysis  was  made  by  Martinelli  et  al  (30)  covering  the  viscous 
flow  region.  These  tests  were  conducted  in  tubes  of  2  to  50  feet  in 

*  Figures  in  parentheses  refer  to  the  bibliography  at  the  end.  (Page  86) 
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length  ranging  in  size  from  capillary  tubing  to  1  inch  pipe.  A  year 
later  in  19^+7 >  Bergelin  and  Gazley  (6)  at  the  University  of  Delaware 
investigated  the  effect  of  the  hydraulic  gradient  in  stratified  flow. 

These  tests  were  carried  out  in  1  inch  pyrex  tube,  the  pressure  taps 
being  10.^8  feet  apart. 

Lockhart  and  Martinelli  (28)  presented  a  final  and  revised 
correlation  in  1949,  covering  all  the  flow  mechanisms.  This  correla¬ 
tion  was  based  not  only  on  the  data  from  their  previous  work  but  also 
one  some  new  data  covering  a  wider  range  of  flow  rates  and  liquid 
properties.  Air  was  the  only  gas  used  though  eight  different  liquids 
including  water,  benzene,  kerosene,  etc.  were  employed.  Martinelli  and 
coworkers  (28,  29,  30 )  basing  their  correlation  upon  an  analogy  to 
single  phase  flow,  suggested  the  possibility  of  the  following  mechanisms 
existing  in  two -phase  flow: 

i)  Both  the  liquid  and  gas  in  turbulent  motion 

ii)  The  liquid  in  viscous  and  gas  in  turbulent  motion 

iii)  Both  the  liquid  and  gas  in  viscous  motion 

iv)  The  liquid  in  turbulent  and  gas  in  viscous  motion 

The  transitional  Reynolds  number  from  viscous  to  turbulent 

flow,  calculated  as  if  each  phase  was  flowing  alone,  was  arbitrarily 
chosen  as  2000  for  the  gas  and  1000  for  the  liquid  phase.  The  correla¬ 
tion  is  also  based  on  the  following  two  basic  assumptions  which  are 
quoted: 

i)  "that  the  static  pressure  drop  for  the  liquid  phase 

must  equal  the  static  pressure  drop  for  the  gas  phase, 
regardless  of  the  details  of  flow  as  along  as  appreciable 
static  pressure  differences  do  not  exist  along  any  pipe 
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ii)  that  the  volume  occupied  by  the  liquid  plus  the  volume 
occupied  by  the  gas  at  any  instant  must  equal  the 
total  volume  of  the  pipe". 

Gazley  (6)  in  his  experiments  on  stratified  flow,  found  that 
the  profile  of  the  liquid  surface  varied  along  the  tube,  thereby  showing 
that  the  pressure  drops  experienced  by  each  phase  were  not  the  same 
when  the  flow  was  stratified.  Gazley  measured  the  pressure  drops  in  the 
gas  phase.  Jenkins  (22)  who  measured  the  pressure  drops  in  the  liquid 
phase  found  that  these  were  as  much  as  100$  higher  than  those  in  the 
gas  phase  for  stratified  flow.  This  was  due  to  the  hydrostatic  head 
difference  in  the  liquid. 

There  is  a  general  agreement  among  various  investigators 
(3,  2,  21,  22,  25)  that  six  or  seven  patterns  of  flow  can  exist  depending 
upon  the  flow  rates  and  a  given  set  of  fluids.  Thus  Martinelli's  distinc¬ 
tion  of  viscous  and  turbulent  flow  in  each  of  the  phases  does  not  in 
itself  describe  the  flow  pattern.  The  flow  patterns  that  are  developed 
when  increasing  amounts  of  gas  are  introduced  into  a  pipe  flowing  full 
of  liquid  are : 

i)  Bubble  Flow:  Flow  in  which  bubbles  of  gas  move  along 
the  upper  part  of  the  pipe  at  approximately  the  same 
velocity  as  the  liquid. 

ii)  Plug  Flow:  in  which  alternate  plugs  of  liquid  and  gas 
move  along  the  upper  part  of  the  pipe.  The  plugs  of  gas 
are  in  effect  elongated  bubbles. 

iii)  Stratified  Flow:  in  which  the  liquid  flows  along  the 
bottom  of  the  pipe  and  the  gas  flows  above  it  over  a 
relatively  smooth  gas-liquid  interface. 
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iv)  V/ave  or  Ripple  Flow:  which  is  similar  to  stratified 
flow  except  that  the  gas  moves  at  a  higher  velocity 
causing  wave  formation  at  the  interface. 

v)  Slug  Flow:  Flow  in  which  waves  are  broken  up  by  the 
more  rapidly  moving  gas  to  form  a  frothy  slug  which 
occupies  the  whole  cross-section  of  the  pipe  and  moves 
at  a  much  greater  velocity  than  the  average  liquid 
velocity.  Violent  fluctuations  result  due  to  the  acce¬ 
leration  experienced  by  the  slug. 

vi)  Film,  Semi -annular  or  Annular  Flow :  a  type  in  which  the 
liquid  flows  in  a  film  around  the  inside  wall  of  the 
pipe  and  the  gas  flows  at  a  higher  velocity  as  a  central 
core. 

vii)  Spray  or  Mist-annular  Flow:  a  flow  resembling  annular 
flow,  occurs  at  much  higher  liquid  rates,  in  which  most 
or  all  of  the  liquid  is  entrained  as  a  spray  by  the  gas. 
Schneider  (34)  has  called  this  type  'Cresting’  flow. 

The  Martinelli  correlation  (28)  does  not  cover  slug  flow  because 
of  its  unsteady  nature  and,  as  pointed  out  by  Gazley  (l8)  excludes  strati¬ 
fied  flow  due  to  the  effect  of  the  hydraulic  gradient. 

The  correlation  of  Martinelli  et  al  (28)  involves  the  use  of  an 


experimentally  determined  variable 


o 


lorg 


which  is  defined  as: 


Eqn. 


(D* 


*  Please  refer  to  the  nomenclature  at  the  end.  (Page  84) 
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The  function 


O 


1  or  g 


is  assumed  to  be  dependent  on  one  other 


variable,  X,  which  is  defined  as: 


X  = 


(AP/AL)  • 


nl/2 


(ap/al) 


s 


Eqn.  (2) 


where  the  pressure  gradients  expressed  are  for  the  liquid  and  gas  phase 
calculated  as  if  each  phase  were  flowing  alone  in  the  pipe. 

Four  different  relationships  between  O  and  X  were  obtained 
for  each  of  the  four  mechanisms  possible  depending  upon  whether  the  two 
phases  were  in  viscous  or  turbulent  motion  as  determined  by  their  super¬ 
ficial  Reynolds  numbers.  Another  variable,  R-  ,  defined  as  the  fraction 
of  pipe  volume  occupied  by  the  liquid  phase,  was  measured  and  was  also 
found  to  be  a  function  of  X. 

The  Martinelli  correlation  (28)  represents  their  data  with  an 

accuracy  of  plus  or  minus  thirty  percent.  Later  investigators  (l8,  3^) 

suggested  that  since  it  was  d)  _  that  entered  directly  into  the 

T  1  or  g 

2  2 

calculations  of  two-phase  pressure  drop,  a  plot  of  versus  X  or  ' 

was  more  appropriate  than  that  of  C)  versus  X  as  the  latter  reduced 
the  spread  of  the  data.  Jenkins  (22)  found  an  accuracy  of  plus  or  minus 
40 $  while  comparing  his  data  with  the  Martinelli  correlation.  He  also 
noticed  that  definite  changes  in  slope  occurred  when  his  data  were 
plotted  on  the  O  versus  X  plot,  with  constant  liquid  flow  rates  as 
parameters.  This  he  concluded,  showed  that  (j)  was  not  a  function  of  X 
alone.  His  argument  was  supported  by  the  fact  that  the  changes  in  slope 
coincided  with  the  visually  observed  changes  in  the  flow  pattern  in  all 


cases . 


Some  determinations  of  inter facial  shear  have  been  reported 
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by  Gazley  (17),  who  studied  the  stratified  flow  of  air-water  mixtures 
in  2.065  inch  diameter,  horizontal  tube.  Mechanical  energy  balances 
were  written  for  the  liquid  and  gas  phases  separately.  It  was  assumed 
that  the  wall  friction  factors  for  both  the  gas  and  liquid  could  be 
defined  in  the  usual  fashion  as  a  function  of  the  Reynolds  number 
orovided  the  appropriate  equivalent  diameters  were  used.  The  inter¬ 
facial  shears,  then,  could  be  calculated  from  the  mechanical  energy 
balances  using  the  measurements  of  the  mass  flow  rates,  fluid  temperatures, 
gas  phase  pressure  drop,  liquid  level  and  the  interfacial  gradient. 

The  two  interfacial  shears  are  defined  as: 


where 


q'  (vs  -  vd 

11  =  2  c0 

r  an^ 


and 


Y  J?ig  pg  ( 7g  "  ^l) 
2  Sc 


ig 


are  the  shear  stresses,  lbs.  per  sq.  ft.  at  the  inter - 


il  '  ig 

face  for  the  liquid  and  the  gas  and  f±1  and  fig  are  the  friction  factors 
at  the  interface  for  the  liquid  and  the  gas.  Gazley  (17)  concluded  that 
the  energies  lost  by  the  gas  phase  and  gained  by  the  liquid  at  the 
interface  were  equal  if  the  interface  was  smooth  but  if  it  was  wavy,  more 
energy  was  lost  by  the  gas  phase  than  appeared  as  work  transmitted  to 
the  liquid  phase. 

Van  Wingen  (4o)  compared  the  pressure  drops  for  flow  of  oil 
and  gas  mixtures  in  a  3  inch  flow  line  with  the  Martinelli  correlation  (28). 
He  found  that  at  lower  gas  rates  the  observed  pressure  drops  were  higher 
than  those  predicted  by  the  Martinelli  correlation  but  such  a  difference 
was  less  than  1  psi  per  1000  foot,  which  was  quite  insignificant  from 
the  practical  point  of  view.  At  higher  gas  rates,  however,  the  diver¬ 


gence  was  in  the  direction  of  conservatism  in  that  the  analytically 
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computed  pressure  drops  were  higher  than  the  actually  observed. 

The  work  of  Johnson  and  Abou-Sabe  (23)  on  two-phase  flow  of 
air-water  mixtures  involved  heat  transfer  also.  However,  some  of  their 
data  from  isothermal  tests  in  a  O.87  inch  inside  diameter  pipe  indicated 
that  more  accurate  correlation  was  obtained  if  curves  similar  to  those 
of  Lockhart  and  Martinelli  were  drawn  through  points  having  the  same 
liquid  flow  rates. 

Schueler  (35)  applied  Martinelli 's  correlation  to  liquid- 
liquid  flow.  He  studied  the  isothermal  flow  of  different  mixtures  of 


carbon  tetrachloride  and  water  and  calculated  the  parameters  <P  and  X. 

^  ±  and  O  g  were  found  to  be  53$  lower  than  the  Martinelli  curves. 
(Here  the  subscript  g  refers  to  water  and  1  to  carbon  tetrachloride) 

Levy  (27)  presented  an  analytical  treatment  of  annular  flow 
in  horizontal  pipe  which  confirmed  Martinelli* s  (28)  choice  of  parameters 
He  concluded,  however,  that  for  a  system  in  which  the  ratio  /  |_tg\ 

Ul/ 

exceeded  the  limit  0.005>Pg  >0.001,  additional  parameters  should  be 

KL 

introduced.  His  basic  assumptions  were: 

i)  an  annular  flow  pattern  with  no  effect  of  gravity 

ii)  the  liquid  velocity  at  the  pipe  was  zero  and  the  gas 
and  liquid  velocities  were  equal  at  the  interface  i.e. 
equal  shear  stresses  at  the  gas-liquid  interface. 

iii)  "the  application  of  reduced  Navier-Stokes  equations  for 
laminar  and  the  one -seventh  power  la w  for  trubulent  flow 

The  theoretical  predictions  of  pressure  drop  were  lower  than 
the  empirically  determined  by  Martinelli  by  a  maximum  deviation  of  20$. 
Predicted  values  of  R^j  the  fraction  of  pipe  volume  occupied  by  the 
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liquid  however,  were  considerably  higher  than  the  empirical  values  and 
depended  upon  the  flow  pattern.  He  also  noted  that  the  discrepancies 
between  his  analysis  and  the  experimental  data  might  be  attributed  to 
the  idealized  geometry,  the  assumption  of  equal  shear  stresses  at  the 
gas  liquid  interface,  the  use  of  constant  properties,  the  unsteady 
state  effects  at  the  interface  and  the  effect  of  flow  pattern. 

An  investigation  of  annular  liquid  flow  with  cocurrent  air 
flow  in  horizontal  tubes  was  made  by  Abramson  (l).  Water  containing 
various  quantities  of  detergent  and  water-ethylene  glycol  solutions 
were  injected  radially  into  horizontal  tubes  of  pyrex  or  lucite  and 
visual  and  photographic  observations  made  upon  the  appearance  of  the 
film  surface  as  a  function  of  a  number  of  variables.  A  transition 
region  of  liquid  flow  from  a  relatively  smooth  to  a  definitely  disturbed 
condition  was  observed  for  the  various  liquids  used.  The  limits  of 
this  transition  region  in  terms  of  liquid  flow  rates  were  determined. 

An  optical  device  developed  to  measure  the  film  thicknesses  was  found 
to  be  insensitive  in  the  flow  range  studied  indicating  that  liquid 
film  thicknesses  were  less  than  0.005  inch.  The  velocities  of  the  film 
surface  disturbances  were  measured  by  means  of  a  high  speed  camera. 
Variations  in  air-mass  velocity  between  30  and  108  per  sec."sq.  ft. 
for  4  inch  pyrex  and  lucite  and  2  inch  lucite  tube  resulted  in  appre¬ 
ciable  change  in  the  water  flow  defining  the  transition  region  which  was 
of  the  order  of  0.04-0.07  lbs.  per  sec. -ft.  However,  the  disturbances 
appeared  smaller  and  more  numerous  with  increasing  air  velocity.  The 
effect  of  increasing  viscosity  and  decreasing  surface  tension  was  to 
increase  the  flow  rate  range  over  which  the  transition  occurred.  The 
Reynolds  number  for  water  corresponding  to  the  transition  range  is 


10 


from  200-400.  He  also  applied  the  Universal  Velocity  profile  analysis 
of  Von  Karman  (13)  and  on  this  basis  the  transition  occurred  for  values 
of  Y^  between  12  and  21  corresponding  to  a  region  of  the  buffer  layer 
5  <c  Y*  <  30.  Y+  is  one  of  two  parameters  of  a  relationship  that 
expresses  the  Universal  Velocity  Profile  for  turbulant  flow. 


U 


Px  y 


•*1 


where 


y  =  distance  from  the  pipe  wall,  ft. 

"tT 


'W  c 


1  =  shear  force  at  the  wall  per  unit  area  lbs.  (F)  per  sq.ft, 

w 

Schneider's  (34)  work  for  his  thesis  (M.Sc.  in  Pet.  Eng.)  is 
important  in  the  sense  that  he  was  the  first  person  to  use  a  long  trans¬ 
parent  test  section  for  his  studies  which  covered  mainly  wave  or  ripple, 
film  or  semi-annular  and  annular  flow.  Twenty  foot  sections  of  'Kraloy' 
plastic  tubing  of  1,  1  l/2  and  2  inch  diameter,  with  an  average  overall 
length  of  100  feet  and  an  average  test  section  length  of  70  feet  were 
employed.  Natural  gas  in  combination  with  water,  kerosene  and  SAE  40 
lubricating  oil  were  used.  His  correlation  is  a  plot  on  logarithmic 

co-ordinates  of  a  pseudo  friction  factor  f,  and  the  ratio  G,  u,  . 

tpg  1  1 

G  '  u 

g 

f.  is  defined  by  the  equation: 
tpg 


tpg  = 


(AP/AL)tp  gc  D  p( 

"2 - 


Eqn.  (3) 


2  G 


g 


This  correlation  represents  his  data  with  an  accuracy  of  +32$ 
to  -24$.  It,  however,  applies  to  those  flow  types  which  exhibit  steady 
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state  conditions  and  where  evaporation  and  kinetic  energy  losses  are 
small.  It  does  not  include  slug  and  stratified  flow. 

Alves  (2)  in  his  work  with  air-water  mixtures  in  a  1  inch 
pipe  line  contactor  found  fairly  good  agreement  with  the  Lockhart  and 
Martinelli  (28)  correlation.  In  general  the  experimental  values  were 
higher  than  the  predicted  for  annular  flow  and  lower  than  the  predicted 
for  the  slug  and  wavy  flow. 

Baker  (3,  b)  presented  data  collected  at  Magnolia  Petroleum 
Company  field  research  laboratory,  for  the  flow  of  gas -oil  mixtures, 
through  b  to  10  inch  diameter  lines.  The  length  of  the  lines  varied 
from  2  to  8  miles.  He  found  that  the  Martinelli  (28)  correlation  was 
inadequate  for  the  larger  diameter  lines  and  concluded  that  the  flow 
pattern  had  a  marked  effect  on  the  pressure  drop.  Separate  equations 
for  different  flow  types  were  proposed  including  one  for  the  case  of 
a  hilly  terrain.  These  equations  are  modified  forms  of  the  Martinelli 
equations  involving  O  ,  X  ,  and  ,  the  mass  velocity  of  the  liquid 
phase  lbs.  per  hr.  -  sq.  ft.  of  total  pipe  cross-section. 

Fried  (l6)  reported  that  the  parameter  ^  ^tt  f°r  non" 
isothermal  flow  of  air-water  mixtures  was  about  equal  to  the  Martinelli 
(28)  curve  provided  the  pressure  drop  measurements  were  corrected  for 
momentum  pressure  drop.  In  making  the  momentum  corrections,  Fried 
considered  the  mixture  flow  at  the  inlet  and  exit  as  being  homogeneous. 
The  data  of  Fried  were  obtained  in  a  0-737  inch  inside  diameter  pipe 
15.7  ft.  long,  at  pressures  up  to  40  psig. 

White  and  Huntington  (4l)  presented  an  empirical  correlation 
as  a  result  of  experiments  on  three  liquids  and  two  gases  in  1,  1  l/2 
and  2  inch  diameter  lines,  with  an  overall  length  of  100  ft.  and  test 
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section  length  of  70  ft.  The  correlation  represents  their  data  with 
a  maximum  deviation  of  +  or  -17$.  It  is  a  plot  on  logarithmic  co¬ 


ordinates  of  the  two  variables  e  and  R  defined  as: 

e 


and 


2gc  *Ptp  P1 


W 


3.6 


1 


W1  +Wg 


R  = 


W 


a. 8 


pg 

pr 


,0.9 


0.1 


1 


0.1 


Eqn.  (4) 


Eqn.  (5) 


This  correlation,  however,  is  limited  to  low  gas  phase  densities,  and 
to  liquids  with  viscosities  of  less  than  120  centipoises.  It  also  does 
not  include  slug  and  stratified  flow. 

Another  correlation  that  has  found  favour  with  later  investi¬ 
gators  is  the  one  presented  by  Chenoweth  and  Martin  (ll).  Their  data 
werecollected  using  air-water  mixtures  in  35  ft.  test  sections  of  1  l/2 
inch  and  3  inch  diameter  galvanized  steel  pipes  at  pressures  from  atmos¬ 
pheric  to  100  psia.  They  found  that  the  Martinelli  correlation  (28) 
was  reasonably  suitable  for  the  data  collected  with  the  1  1/2  inch  line 
but  predicted  high  pressure  drops  for  the  3  inch  line.  In  some  cases 
the  predicted  pressure  drop  was  2.5  times  as  great  as  the  actually 
observed.  Their  correlation  is  a  plot  on  logarithmic  co-ordinates  of 
the  superficial  liquid  volume  fraction,  L.V.F. ,  as  abscissa,  the  ratio 

^^tp  >  as  ordinate  with  the  ratio  APr.  as  the  parameter.  The 

ap2*  AP-j  * 

parameter  is  a  ratio  of  a  fictitious  all  gas  pressure  drop  calculated 
from  total  flow  rate  but  using  the  properties  of  the  gas  to  a  pressure 
drop  calculated  in  the  similar  way  but  using  the  properties  of  the  liquid. 
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The  superficial  liquid  volume  fraction  is  calculated  from  the  flow 
rates  and  the  densities  of  the  two  phases  as: 


l.v.f.  =  w.y.p.1 _ 

Wl/Pl  +  Wg/Pg 


Qi  +  Q 


=  c 


L 


g 


It  is  claimed  that  this  correlation  can  he  used  for  any  liquid-vapor 

mixture  as  long  as  the  Reynolds  number  calculated  from  the  total  flow 

rates  and  the  properties  of  the  liquid  is  greater  than  2000.  However 

it  does  not  specifically  include  the  concept  of  flow  pattern. 

Bertuzzi  et  al  (7)  presented  another  correlation  developed 

from  data  then  available  in  literature.  A  friction  factor  f^  is 

tpm 

defined  as: 


'tpm 


(AP)tp  Bc  D  Pm 

2  L  G  2 
m 


Eqn.  (6) 


In  the  correlation,  f ,  is  plotted  as  a  function  of  the  product 

a  13 

(Nr,.),-.  ( Np,e )  on  logarithmic  co-ordinates. 


The  exponents  a  and  b  are  defined  as  a 


R 


M 


and  b  = 


%fl 


o.m 


M 


Where  R^  is  the  gas-liquid  mass  ratio. 

The  selection  of  the  exponents  a  and  b  was  arbitrary  except  that  they 

met  requirements  at  the  limiting  conditions  of  all  liquid  and  all  gas. 

The  plot  was  divided  into  four  groups  with  a  range  of  values  of  gas -liquid 

mass  ratio  as  parameter  i.e.  R  values  of  0-0.2,  0. 2-0.4,  0.4-0. 6  and 

M 

0. 6-1.0.  Furthermore,  the  individual  correlation  for  a  particular  range 
of  Rj.  values  consisted  of  two  distinct  curves  which  gave  two  values  of 
f  for  a  single  value  of  U‘^e)0.  (%e)i  •  ®ne  curve  stated  to 

apply  to  stratified,  wave  and  semi -annular  flow  and  the  other  to  other 
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types.  The  transition  from  one  curve  to  the  other  appeared  to  "be 
associated  with  liquid  Reynolds  numbers  in  the  region  of  500  to  10,000. 

Reid  et  al  (32)  conducted  tests  in  4  inch  and  6  inch  commercial 
steel  pipes  with  an  overall  length  of  76  ft.  and  with  the  pressure  taps 
56  ft.  apart  to  determine  the  validity  of  Chenoweth  and  Martin  (ll)  and 
Lockhart-Martinelli  (28)  correlations  for  low  pressure  air-water  flow 
in  large  diameter  pipes.  It  is  noted  that  the  data  from  the  4  inch  line 
were  well  represented  by  the  Lockhart-Martinelli  correlation  while  the  6 
inch  data  fell  about  20 $  below  but  on  a  smooth  parallel  curve.  Comparison 
of  experimental  data  with  Chenoweth  and  Martin  (ll)  correlation  showed 
excellent  agreement  and  no  noticeable  effects  of  pipe  diameter.  They 
concluded  that  above  a  superficial  liquid  volume  fraction  of  0.1  i.e.  10 $ 
to  100$,  two-phase  pressure  drop  may  be  estimated  from  the  equation: 

APtp  =  AP*  (LVF)"°‘U'J  Eqn.  (7) 


since  at  high  values  of  superficial  LVF  ( >  0.1)  all  parameteric  curves 
converge  to  a  single  line. 

Chisholm  and  Laird  (12)  presented  two-phase  air-water  flow  data 
in  one  inch  diameter  rough  tubes.  Besides  a  smooth  brass  tube,  five  other 
tubes  were  employed  with  the  £/D  value  ranging  from  0.0025  to  0.068. 

The  test  section  was  8  feet  long  with  the  pressure  taps  located  at  2 
feet  intervals.  The  tubes  were  in  the  horizontal  plane  and  the  pressures 
were  close  to  atmospheric  throughout.  According  to  these  investigators 
the  general  form  of  Lockhart  and  Martinelli  (28)  parameter: 
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n,  being  the  power  of  the  Reynolds  number  in  the  friction  factor 
relationship  f’  =  C'  N^e  n  Eqn.  (9) 

was  not  very  appropriate.  They  argued  that  though  for  smooth  tubes, 
for  which  the  Blasius  Equation  gave  n  =  0.25;  Eqn.  8  was  reduced  to: 


Eqn.  (l0) 


and  provided  a  satisfactory  parameter  for  the  correlation  of  two -phase 
flow  data,  it  did  not  provide  a  satisfactory  parameter  for  the  corre¬ 
lation  of  two -phase  flow  data  in  rough  tubes  for  which  the  value  of  n 
approaches  zero.  Moreover,  they  found  that  the  most  satisfactory  corre¬ 
lation  for  rough  tubes  data  was  obtained  using  X  calculated,  regardless 
of  the  value  of  n,  by  Eqn.  (l0).  Hence  they  redefined  the  parameter  as: 


875  ^'°'125  —°-5 
M- 


Eqn.  (ll) 


In  the  special  case  of  n  =  0.25,  X  is  identical  to  X,  the  Martinelli 
parameter. 

AP 

They  recommended  logarithmic  plots  of  tp  -1,  as  a  function 


AP 


1 


AP, 


of  parameter  X  in  preference  to  the  plot  of  tp  as  a  function  of 

APX 

the  parameter  X  as  recommended  by  Lockhart  and  Martinelli.  For  rough  tubes 
the  following  equation  is  recommended: 


AP 


tp 


1  + 


APn 


m 


Eqn.  (12) 


Where  C  and  m  are  constants.  Approximate  correlations  for  the  estimation 


of  C  and  m  were  given  in  the  form  of  arithmetic  plots,  wherein  C  was 
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plotted  as  a  function  of  J(%e),  ^RT  and  '  ni '  as  a  function  of  '"RT  , 

xs  xs 

where  (%e)l  =  Reynolds  number  for  the  liquid  flowing  alone 

dimensionless. 

fR«p  =  friction  factor  for  rough  tubes,  dimensionless. 

fs  =  friction  factor  for  smooth  tubes,  dimensionless. 

This  correlation  gave  pressure  drop  predictions  within  plus  or  minus  15$ 
of  the  experimental  values  for  the  majority  of  their  tests. 

They  presented  another  equation  that  correlated  the  two -phase 
pressure  drops  and  'saturations'.  Saturation  was  defined  as  the  fraction 
of  the  pipe  volume  occupied  by  each  phase. 


AP. 


1.6  f  L  (V-J  )  P 
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tp 
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Eqn.  (13) 


Where  Vj  =  Actual  lineal  velocity  of  the  liquid  during  two-phase 

flow  ft./ sec.  and  is  equal  to  _1  wheieV-,  is  the  superficial  velocity 

R1 


of  the  liquid  phase  based  on  the  total  pipe  cross  section  and  R^  is  the 
liquid  saturation,  the  volume  fraction  of  pipe  occupied  by  the  liquid. 

The  correlation  in  the  form  of  equation  (l3),  however,  is  limited  to 
commercial  tubes  with  ^/D  values  up  to  0.0025.  No  satisfactory  correlation 
of  this  form  could  be  obtained  with  the  rougher  tubes  since  in  these  cases 
the  data  no  longer  fell  on  a  single  line.  They  also  developed  equations 
for  the  prediction  of  liquid  saturation  having  the  general  form: 
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a,  b,  c  and  d  are  all  constants. 
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Besides  the  North  American  continent,  considerable  work  is 
being  done  in  this  field  in  other  parts  of  the  world,  particularly 
in  Europe  and  the  U.S.S.R. 

Kosterin  (24)  showed  that  the  friction  factor  (hydraulic 

resistance)  in  the  flow  of  gas -liquid  mixtures  through  horizontal 

pipes,  was  a  variable  quantity  depending  on  the  concentration  of  the 

gas  phase  (vapor  content)  the  Reynolds  number  calculated  for  the  liquid 

phase,  the  diameter  and  roughness  of  the  pipe,  the  pressure  at  small 

Reynolds  numbers  for  the  liquid.  Comparison  with  experiments  performed 

on  water-vapor  mixtures  at  high  pressures  showed  the  analogy  with  the 

hydrodynamic  relationship  obtained  from  experiments  on  air -water  mixtures 

The  same  investigator  (25)  reported  the  results  of  experiments 

with  air-water  mixtures  in  1,  2,  3  and  4  inch  diameter  tubes  with  49 

foot  test  sections.  He  found  that  the  friction  factor  was  a  function  of 

the  volumetric  flow  concentration  of  the  gas  and  the  Froude  number. 

He  defined  the  friction  factor  f  for  two -phase  flow  as: 

tpm 


AP 


tp 


1  D  p  A 
'c  m 


2  AL  W 


Eqn.  (l4) 


m 


Where  p  =  density  of  the  mixture,  lbs 
m  - 


and  is  calculated  as : 


cu.  ft. 


m 


p  C  +  (l-C  )  p  where 
g  S  g  1 


Eqn.  (15) 


C  =  Volume  flow  concentration  of  the  gas  and  is  given  by 
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Q  +  Q. 


where  the  Qs  are  the  volumetric  flow  rate  of  each  component,  cu.  ft. 
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The  coefficient 
of  the  two-phase  friction  factor  f 


ffM  ^  ,  representing  the  deviation 


tpra  ^Be'tpl 


tpm 


from  the  law  of  resistance  w as 


defined  as  the  ratio  of  f  to  the  friction  factor  of  a  homogeneous 

tpm 

liquid  taken  for  the  Reynolds  number  (®pe)tpl  =  ^ 


jt  D  n. 


The  Froude  number  w as  defined  as : 


W 


Fr_ 


J1L 


m 


.2  2 
g  A  Dp 
c  m 


Eqn.  (l6) 


An  arithmetic  plot  of  the  coefficient  °  against  C  with  Fr  as  para- 

g  m 

meter  is  given.  The  range  of  Fr^  is  from  0.1  to  1,000. 

Kosterin  and  Rubanovich  (26)  carried  out  tests  to  determine  the 
effects  of  surface  tension  on  pressure  losses  in  two-phase  flow.  It 
was  reported  that  a  decrease  in  the  surface  tension  of  the  liquid  phase 
did  not  exert  an  influence  on  the  magnitude  of  pressure  losses  under 
conditions  of  unstable  foam  formation  over  a  considerable  range 
(0.164  to  11.481  ft.  per  sec.)  of  velocity  of  the  liquid  phase.  A 
decrease  in  surface  tension  however,  increased  the  degree  of  dispersion 
of  gas  in  the  liquid  in  certain  types  of  flow.  The  formation  of  an 
intense  and  stable  foam  in  pipe  could  lead  to  a  considerable  increase 
in  the  'hydraulic  resistance'  during  the  flow  of  the  gas-liquid  mixtures. 

Hoogendoorn  (21)  employed  both  smooth  pipes  with  the  inside 
diameter  ranging  from  0.945  inch  to  5*51  inch  and  rough  pipe  of  1*97  inch 
inside  diameter  with  the  c /D  value  ranging  from  0.0012  -  0.03,  in  his  experi¬ 
ments.  Water,  spindle  oil,  and  gas  oil  in  combination  with  air  were 
used  in  a  horizontal  82  foot  long  test  section  of  transparent  perspex 
pipe.  He  found  that  the  Lockhart-Martinelli  (28)  correlation  cannot 
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be  used  for  the  following  three  flow  patterns: 

a)  Plug,  slug  and  froth  flow,  if  gas  densities  differ  from 
that  of  air  at  atmospheric  pressure. 

b)  Stratified  and  wave  flow. 

c)  Mist -annular  flow. 

Separate  correlations  for  each  of  these  patterns  were  suggested. 
For  plug,  slug  and  froth  flow,  an  experimental  correlation  which  included 
the  influence  of  gas  density  and  pipe  roughness  gave: 
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Eqn.  (IT) 


Eqn.  (l8) 


/G  \ 

_g  <  0.05  and  (%e)^p3_  >  3000,  where 

\  m  ' 


the  limiting  conditions  being 
* 

AP.  =  Pressure  drop  over  distane  AL,  calculated  for  the  total 

mass  (W-  +  Wn  )  of  gas  and  liquid,  flowing  in  the  liquid  phase. 

(N„  V  1  =  GmD  where  u  is  the  viscosity  of  the  liquid  in  — - 

tfe  tpl  - —  1  ft.  sec, 

^1 


For  wave  flow  a  new  correlation  gave: 
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Eqn.  (19) 


Where  C  is  a  constant,  whose  value  depends  mainly  on  pipe 
diameter  and  pipe  roughness.  Some  values  of  C  were  given.  For  mist- 
annular  flow,  he  found  that  the  pressure  drop  was  independent  of  the 
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liquid  mass  velocity  if  it  was  greater  than  6.l4  lb . _  j 

sq.  ft.  -  sec. 


’■Okg. 


sq.  m. 


-  sec 
/ 


A  new  correlation  based  on  data  from  air-gas  oil  runs  gave: 
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Eqn.  (20) 
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The  following  correlation  was  given  between  the  slip  velocity ,  V  ,  and 

S  jL 


the  gas  hold-up,  R  : 
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Rg  =  0.6  |  R  V  .  1°'°5 
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Eqn.  (21) 


Eqn.  (22) 


Eqn.  (21)  can  be  written  as 
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g  -1 

Eqn.  (23) 


This  correlation,  however,  is  based  only  on  the  air-oil  flow  data. 

Widell  (4l)  conducted  tests  with  air-water  mixtures  in  0.984 
inch,  1.4l6  inch  and  3*36  inch  inside  diameter  brass  pipes  and  with  steam 
and  water  in  1.008  inch  copper  pipe.  The  distance  between  the  pressure 
taps  was  very  small  being  3*395  ft.,  4.92  ft.,  5*9  ft.  and  3*34  ft. 
respectively.  The  Martinelli  correlation  (28)  was  found  to  be  unsuitable 
for  his  data  especially  at  lower  water  velocities.  No  correlation  in 
the  real  sense  was  presented  although  logarithmic  plots  of  vs 
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G  (G^  +  Gc.)  were  given  for  the  0.984  inch  and  1.4l6  inch  diameter  pipes, 
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Baxendell  (5)  has  derived  an  empirical  relationship  linking 
density,  velocity  and  pipe  diameter  to  predict  pressure  losses  in 
two -phase  pipe  line  flow.  The  final  equation  was  the  result  of  an 
analysis  of  each  term  involved  in  the  integrated  form  of  the  mechanical 
energy  balance  for  flow  of  a  compressible  fluid  in  a  circular  pipe  of 
constant  diameter  as  applied  to  two-phase  flow,  particularly  oil  and  gas. 
It  is  assumed  that  the  density  of  the  mixture  might  be  expressed  in 
terms  of  pressure  through  the  gas  law. 


The  final  equation  is: 
/2 


-2 


2  2  9.187  F  Q’  L’ 

pap  =  px  -p2  -  - -2-5 - 

m  K  D 


Eqn.  (24) 


where  Q' 
L' 
D 
m 


=  rate  of  flow  of  tank  oil 


K 

Y 

X' 


Pg 


cu.m. 

day 


length  of  pipe  line  in  kilometers 
diameter  of  pipe,  ft. 

slope  of  the  density-pressure  curves  at  various  gas¬ 
oil  ratios,  sq-* 

cu.ft. 


(24)  (3600) 


where 


(Y  +  X'  Pg) 

mass  of  one  cu.m,  tank  oil  +  soln.  gas,  lbs. 
excess  gas-oil  ratio  cu.m. /cu.m. 


=  density  of  gas  at  standard  conditions, 


=  friction  factor 


lbs, 

cu.m. 


Another  relationship  that  is  necessary  for  the  calculations  is: 
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L  or  K  = 
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WT  is  the  total  mass  rate  of  flow  -  sum  of  the  mass  rate  of 


flow  of  tank  oil  plus  the  mass  rate  of  flow  of  gas, 


Since  the  terms,  'm'  and  'K'  are  characteristics  of  a  particular 
field,  it  is  necessary  that  these  be  evaluated  before  equation  (24)  can 
be  used  to  calculate  the  pressure  drop.  For  field  calculations,  the 
first  step  is  to  construct  from  the  PVT  data  of  a  representative  well 
a  set  of  pres sure -density  curves  to  cover  the  range  of  gas-oil  ratios 
encountered  in  the  area.  From  these  curves,  plots  of  gas-oil  ratios 
against  m  and  K  are  drawn.  A  plot  of  friction  factor  F  against  DVP, 

0 ' 

which  in  terms  of  field  units  can  be  written  as  equivalent  to  1.2733  Spr 

jSD 

or  1.2733  Wp  ,  is  also  necessary.  This  plot  is  constructed  from  the 
D 

field  data.  At  the  high  rates  of  flow  generally  encountered  in  pipe 
lines,  the  degree  of  turbulence  is  of  such  magnitude  that  of  the  total 
energy  loss,  the  portion  due  to  viscous  shear  is  negligible.  Consequently, 
in  the  turbulent  zone,  a  plot  of  friction  factor  against  DVP  can  be 
expected  to  give  a  definite  trend  equivalent  to  the  curves  of  Reynolds 
number  versus  friction  factor  for  single  phase  fluid  flow. 


Freeman  (15)  has  also  developed  an  empirical  relationship 


for  two -phase  flow,  from  the  data  obtained  at  the  Kirkuk  oil  field  in 
Iraq. 


Eqn.  (25) 


where 


P  =  Pressure,  in  ■  ^  abs. 

in. 


L  =  Length  1000  ft 


Q  =  Flow  rate 


1000  bbl. 
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Diameter  inches. 
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Teletov  (39)  defined  a  friction  factor  similar  to  the  one 
employed  by  Kosterin  (25)  i.e. 


=  t  ^ 


) 


Re ' tpm 


tpm 

where  the  function  CL>  representing  the  deviation  of  the  two-phase 
friction  factor,  f-tom»  from  the  law  of  resistance  was  defined  as  the 
ratio  of  ftpm  the  friction  factor  of  homogeneous  liquid  taken  for 

z  .  k  u 

the  Reynolds  number  C^Re  )tpm  =  ' >  where 


it. 


ll  =  Effective  viscosity  of  the  mixture,  — — — 

ft. sec. 

The  effective  viscosity  of  the  mixture  was  defined  using  the  corresponding 
kinematic  viscosities  as  follows: 


2/, 


_g 

^g 


+  ci 


g 


^1 


m  yg  UL 

Two  arithmetic  plots  of  and  with  the  viscosity  ratio 
as  the  parameter,  for  viscous  and  turbulent  flow  were  given. 


Eqn.  (2 6) 


A  considerable  amount  of  work  has  been  done  at  the  University 
of  Alberta,  on  two -phase  flow  both  in  vertical  and  horizontal  pipes. 
Studies  of  vertical  gas-liquid  flow  at  this  university,  by  Govier  and 
his  associates  have  led  to  a  separation  of  the  hydrostatic  head  component 
from  the  irreversibility  component.  This  was  accomplished  by  Radford  (31) 
who  investigated  the  effect  of  air  and  water  flow  rates  on  the  flow 
pattern.  By  the  application  of  a  mechanical  energy  balance  to  each  of 
the  two  flowing  fluids,  he  was  able  to  develop  a  general  equation 
which  permitted  the  separation  of  the  pressure  drop  into  its  two 
components.  The  equation  is  as  follows: 
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V1  ,  1+JVj 


+ 


AX 


1  +  R 


1  +  R 


v 


AF  ) 
AX/ 


Eqn.  (27) 


/a F  ) 

where!  /  =  Sum  of  the  irreversibility  effects  attributable  to 

both  the  gas  and  liquid  phases  in  height  AX,  :.(i. 

I  X>! 

1  +  Rw 

It  can  be  seen  from  the  equation  (27)  that  the  term  - — — —4  represents 


1  +  R. 


v 


the  hydrostatic  head  component  and  the  term 


AF 


1  +  Rv  \AX 


represents 


the  irreversibility  component.  Moreover,  in  the  limiting  conditions  of 
al  1  liquid  flow,  equation  (27)  reduces  to: 


-v. 


AP 

AX 


—  1  + 


AF 

sz 


Eqn.  (28) 


where 


AF  _  2  f  V4 

AX 


gc° 


since  then  R^  =  R^  =  0  and  where  f  is  the  common 


single  phase  friction  factor  and  V,  the  fluid  velocity,  in  feet  per  sec. 

Dunn  (l4)  continued  this  work  and  was  able  to  determine  the 
effect  of  the  air  and  water  flow  rates  on  the  pressure  drop,  flow  pattern 
and  hold-up.  He  proposed  a  correlation  based  on  the  flowing  gas  phase, 
using  a  superficial  friction  factor  and  Reynolds  number.  Results  of  the 
work  of  Radford  and  Dunn  were  published  by  Govier,  Radford  and  Dunn  (l9 ) • 
Their  correlation  is  a  logarithmic  plot  of  the  friction  factor  f  ^  ^  > 
based  on  the  liquid  properties  and  superficial  liquid  velocity  V-^,  with 
the  gas-liquid  input  volume  ratio  Rv  as  the  parameter.  The  friction 
factor  is  defined  by  the  following  equation: 
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Eqn.  (29) 
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Another  correlation  represents  the  hold-up  ratio  IL^  as  a  function  of 
the  superficial  liquid  velocity  VL,  with  Ry  again  being  the  parameter. 
The  hold-up  ratio  HR,  is  defined  as  the  ratio  of  (a)  the  gas -liquid 
volume  ratio  in  the  supply  mixture  to  (b)  the  gas-liquid  volume  ratio 
in  the  flow  section,  all  volumes  calculated  at  average  flow  conditions. 

They  also  observed  that  a  plot  on  arithmetic  co-ordinates 

V  Ap 

of  unit  pressure  drop,  1  tp  ,  and  gas-liquid  volume  ratio  R_  at  a 

“ax”  v 

constant  liquid  rate,  could,  at  low  liquid  velocities  be  divided  into 
four  sections,  which  they  called  the  "Pressure  drop  Regimes".  The 
following  flow  patters  were  observed  both  visually  and  photographically. 

l)  Bubble  2)  Slug  3)  Froth  4)  Ripple  5)  Film  or 
annular  and  6)  Mist. 

A  closer  look  at  a  typical  pressure  drop  relationship  at 
low  liquid  rates  shows  that  the  bubble  flow  pattern  extends  to  the 
first  half  of  Regime  I,  the  slug  pattern  occurs  during  the  second  half 
of  Regime  I  and  extends  midway  into  Regime  II.  The  second  half  of 
Regime  II  and  most  of  Regime  III  is  occupied  by  froth  flow.  Ripple, 
film  or  annular  and  mist  flow  occur  in  the  last  Regime.  The  above  flow 
patterns,  however,  do  not  occur  for  all  water  rates.  As  the  water  rate 
is  increased  Regime  II  and  III  become  progressively  narrower  and  finally 
disappear. 

In  1957,  Short  (36)  expanded  the  experimental  facilities  used 
by  Radford  (31)  and  Dunn  (l4)  and  was  able  to  determine  the  effect  of 
tube  diameter  on  the  pressure  drop,  flow  patterns  and  hold-up.  Data 
collected  by  Short  (36)  together  with  the  applicable  data  from  the  work 
of  Radford  (31)  and  Dunn  (l4)  was  published  by  Govier  and  Short  (20). 
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New  correlations  are  presented  which  enable  the  prediction  for  any 
diameter  between  0.5  and  3*0  inch,  of  the  pressure  drop,  the  hold-up 
and  the  flow  pattern,  as  a  function  of  the  diameter  and  the  gas  and 
liquid  rates  for  the  upward  vertical  flow  of  air-water  mixtures  at  an 
average  pressure  and  temperature  of  36. 0  psia  and  70-85°F  respectively. 

The  friction  factor  f '  ,  is  correlated  as  a  function  of  a 

tpl 

reduced  form  of  the  superficial  Reynolds  number  for  the  liquid  phase 
DV^,  with  as  the  parameter. 

Three  different  correlations  for  the  hold-up  ratio,  Hp,  one 
each  for  Regimes  I  and  II,  III  and  IV  are  presented.  The  hold-up  ratio, 

Hp  ,  and  the  superficial  water  velocity,  V  ,  are  the  ordinates  and  abscissae 
respectively  in  each  case,  the  parameter  being  RD^  3 }  and  R^D  ‘*"/  3 
respectively. 

Further  studies  of  two-phase  vertical  flow  were  made  by 
Sullivan  (38)*  The  two-phase  systems  studied  were  composed  of  water  as 
the  reference  phase  with  air  or  oil  as  the  second  phase.  His  work  was 
an  investigation  of  the  effect  of  the  second  phase  (non -water)  density 
on  the  pressure  drop,  flow  patterns  am  the  hold-up.  Three  second 
phase  densities  were  studied  by  varying  the  average  pressure  in  the 
air-water  system  and  a  fourth  was  obtained  with  oil  as  the  second  phase. 

The  friction  factor,  fp,  correlation  is  a  logarithmic  plot  of 
fp  and  DVp  with  i)3*85ysl.l(2^-l^)0*33  as  the  parameter,  where 
fp  is  the  friction  factor,  dimensionless,  calculated  from  the  equation 
(29)  developed  by  Govier,  Radford  and  Dunn  (19).  The  subscript  'R' 
refers  to  the  reference  phase,  defined  as  the  phase  having  the  lower 
kinematic  viscosity  (water  here)  and  ’s’  to  the  second  phase  (air  or  oil). 
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The  correlation  for  the  hold-up  ratio  H_,  is  a  semi-log  plot 
of  !I^  VR°-0022  ^^R  ps^  and  pp  -  pn  with  J  as  the  parameter. 

Data  collected  by  Sullivan  with  the  air-water  system  together 
with  applicable  data  from  previous  investigators  in  this  field  are  being 
published  by  Brown,  Sullivan  and  Govier  (9). 

Studies  of  two-phase  flow  in  horizontal  pipes  at  the  University 
of  Alberta  and  the  Research  Council  of  Alberta,  have  been  centered  around 
only  liquid-liquid  flow.  Russell  (33)  studied  the  flow  characteristics 
of  a  white  mineral  oil  'Kremol  70'  (having  a  specific  gravity  of  0.83*1- 
and  viscosity  of  22  c.p.  at  20°C)  in  the  presence  of  water  in  a  smooth 
one  inch  pipe.  The  pressure  drop  data  are  presented  as  a  modified 
Fanning  friction  factor,  f-tri,  (based  upon  water  properties)  versus  a 
superficial  velocity  of  water,  V^,  with  the  input  volume  ratio  of  oil  to 
water  as  the  parameter.  The  friction  factor  f^,_^  is  defined  as: 


'tpl 


AptP  gc  D 

2  L  Vp2  Px 


Eqn.  (30) 


(The  subscript  ' 1 '  refers  to  the  water  phase) 

Three  flow  patterns  were  observed  which  were  designated  as 
bubble  flow,  stratified  flow  and  mixed  flow.  Hold-up  measurements 
indicated  that  in  the  laminar  region  of  flow,  the  hold-up  ratio,  H^, 
was  a  function  of  the  input  volume  ratio  and  the  viscosities  of  the 
two  liquids,  while  in  the  turbulent  region,  was  also  a  function  of 
the  superficial  water  velocity. 

Charles  (10)  studied  the  flow  of  three  oils  of  viscosities 
6.29,  16.8  and  65.O  cp  in  the  presence  of  water  in  a  horizontal  pipe 
line  of  1.04  inch  internal  diameter.  The  densities  of  the  oils  were 
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equalized  to  that  of  water  by  adding  carbon  tetrachloride.  It  was 
noticeable  that  the  oil  viscosity  had  little  or  no  effect  on  the  flow 
pattern,  although  the  most  viscous  oil  showed  dissimilar  behaviour  at 
higher  input  oil-water  ratios.  A  similar  trend  was  observed  with  regard 
to  the  pipe  line  pressure  gradients  especially  at  low  oil-water  ratios. 
At  high  oil-water  ratios,  however,  the  most  viscous  oil  gave  relatively 
high  pressure  gradients  in  those  regions  of  flow  in  which  oil  was  the 
continuous  phase  in  contact  with  the  pipe  wall. 

The  hold-up  relationships  at  constant  superficial  water  velo¬ 
cities,  exhibited  a  maximum  at  input  volume  ratios  corresponding  to 
'oil  in  water  concentric'  flow  regime. 

Theory 

1. 

fluid  in 
the  well 
defined : 


where 

AL  =  distance  across  which  the  pressure  drop  occurs,  ft. 

f  =  single  phase  friction  factor,  dimensionless 

V  =  fluid  velocity  at  average  flow  conditions,  ft.  per  sec. 

p  =  fluid  density  at  average  flow  conditions,  lbs.  per  cu.ft. 

g  =  dimensional  conversion  constant,  lb.  mass  ft . 

lb.  force  sec.2 


Pressure  Drop  and  Friction  Factor: 

The  pressure  drop  and  the  friction  factor  for  the  flow  of  a 
one  phase  through  a  horizontal  pipe  line  are  interrelated  by 
known  Fanning  Equation  through  which  the  friction  factor  is 

AP  2  f  V2  p 

AL  "  “i^D  E(ln-  31) 

AP  =  pressure  drop,  lbs.  per  sq.  ft. 
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pipe  diameter,  ft. 
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Extensive  theoretical  and  experimental  work  has  shown  that  the  friction 
factor  is  a  function  of  the  Reynolds  number,  DVp,  and  the  relative 


roughness,  of  the  pipe  wall. 
D 


The  Fanning  equation  may  be  applied  in  any  of  several  ways 
to  the  flow  of  two  phases  and  serves  to  define  "two-phase  friction 
factors"  as  follows: 

Based  on  the  superficial  velocity  and  density  of  liquid  phase 


APtp  _  g  f  tpl  vl2  0X  _  2  ftpl  °12 

AL  gc  D  gc  D  Pi 


Eqn.  (32) 


Based  on  the  superficial  velocity  and  density  of  gas  phase 


AP. 


tp 


AL 


2  ftpg  vg2  PS 


2  f  G„c 
tpg  S 

D  Pg 


Eqn.  (33) 


Based  on  the  superficial  velocity  and  density  of  the  mixture 


where 


AP 


tp 


AL 


AP 


tp 


AL 


tpl 


f 

tpg 


f 


tpm 


2  f  V  2  p 
tpm  ra  rn 


2  f  G  c 
tpm  m 

gc  D  Pm 
^  m 


Eqn.  (34) 


Two-phase  pressure  drop  over  a  distance  AL  ft.  lb, 


sq.  ft. -ft, 


Two -phase  friction  factor  based  on  the  liquid  phase 
properties,  dimensionless. 


Two-phase  friction  factor  based  on  the  gas  phase 
properties,  dimensionless. 

T  o-phase  friction  factor  based  on  the  properties  of 
the  mixture,  dimensionless. 
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V 


v  , 

nr 


pl' 


V 

P> 
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superficial  velocity  of  the  liquid  at  average  flow  conditions  and 
based  on  the  total  cross-section,  ft.  per  sec. 
superficial  velocity  of  the  gas  at  average  flow  conditions  and 
based  on  the  total  cross-section,  ft.  per  sec. 

superficial  velocity  of  the  mixture  at  average  flow  conditions  and 
based  on  the  total  cross-section,  ft.  per  sec. 

density  of  the  liquid  at  average  flow  conditions,  lb.  per  cu.  ft. 
density  of  the  gas  at  average  flow  conditions,  lb.  per  cu.  ft. 
density  of  the  mixture  at  average  flow  conditions,  lb.  per  cu.  ft. 


Gj_,  mass  velocity  of  the  liquid  phase 


lbs. 


sec.  square  ft. 


based  on  the 


total  cross-section 


V 

G  , 
m 


gc> 


lbs 


mass  velocity  of  the  gas  phase  based  on  total  cross-section  — 

S6C  •  SQL  •  iti 

mass  velocity  of  the  mixture  based  on  total  cross-section  lbs . 


sec.  sq.  ft. 


dimensional  conversion  constant 


lb.  mass 
lb. force 


ft. 


sec . 


D,  diameter  of  the  pipe  ft. 


The  three  two-phase  friction  factors  so  defined  are  mathematically 
interrelated  and  the  choice  among  them  can  only  be  arbitrary. 

Russell  (33)  and  Charles  (10)  have  used  the  friction 
factor,  f^ol,  defined  by  Eqn.  (32)  in  their  work  on  oil-water 
flow  in  horizontal  pipe.  Schneider  (3*0  has  used,  ,  defined 

by  Eqn.  (33)*  Lockhart  and  Martinelli  (28)  have  introduced  a 

p 

parameter,  ^  >  which  is  the  ratio  of  the  two-phase  pressure  drop  to 
the  single  phase  pressure  drop  of  either  the  gas  or  liquid  i.e. 
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AP 


tp  = 


AL 


'  g 


AP 


g  = 


AL 


i 


f 2 

AL 


■  f 

g 


2  ^  2 

2  f  G 
g  g 


ScDPg 


■4; 


2  f  G  2 
1  1 

gc  D  P1 


Eqn.  (35) 


where  f  and  f  are  the  single  phase  Fanning  friction  factors  for  gas 
8  1 

and  liquid  respectively. 

Comparison  of  Eqn.  (35)  with  Eqn.  (32)  and  Eqn.  (33)  shows 


that  f 


tpg 


=  t  and  f  = 

Tg  V  tpi  i 


.  Another  parameter  introduced  by 


Lockhart  and  Martinelli  (28)  is  X  which  is  defined  as  the  ratio  of 
the  single  phase  liquid  pressure  drop  to  the  single  phase  gas  pressure 
drop  i.e. 

2  AD  /  at  ~  _  _  2 


=  ^  =  2  fi  °i  /ee  d  px 

^  2  G  2/g  D  P 

s  g  c  g 


Eqn.  (36) 


X 


g 


1 


1' 


Eqn.  (36) 


Eqn.  (37) 


and  C)  is  plotted  as  a  function  of  X  in  the  general  correlation. 
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Chenoweth  and  Martin  (ll)  have  used  the  ratio  ^tp 


* 


AL 


AP 

AL 


as  the  ordinate  in  their  correlation  where 


AP 

AL 


is  a  fictitious  all 


liquid  pressure  drop  per  unit  length  calculated  for  the  total  mass  flow 
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rate  of  both  liquid  and  gas  but  using  the  properties  of  the  liquid. 
This  again  can  be  written  in  terms  of  a  modified  Fanning  equation  as; 


2  f .  G  2/g  D  p 

tpg  g  '  &c  g 


m 


(AP',Al)1  S  ^  G^/gc  D  px 


2  f  ,  G  J/g  D  p, 

tpl  1  7  C  K1 

2  4  Gm/sc  D  P1 


Eqn .  ( 38 ) 
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Ig  / 
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g 

1 

V  m' 

or 


tpg 


=  C  f 


mi  I  m 

G 


p 


Eqn.  (39) 


where  C  is  the  ratio  of  the  two-phase  pressure  drop  to  the  all  liquid 
m 

t 

pressure  drop  and  f^  is  the  single  phase  friction  factor  calculated 
from  the  total  mass  flow  of  both  the  liquid  and  gas  and  using  the  liquid 
properties . 

Bertuzi  et  al  (7)  and  Kosterin  (25)  have  based  the  two-phase 
friction  factor  on  the  total  mass  velocity  of  the  mixture  and  the  mixture 
density  i.e. 


AP 


tp 


AL 


2  f  G  2 
tpm  m 

gc  D  pm 


Eqn.  (34) 


and  the  mixture  density,  p_  is  defined  as  pT  =  C_  p  +  (l-C  )  pn  where  C 

111  ra  g  g  g  1  g 

is  the  volume  fraction  of  the  gas  in  the  mixture.  Comparing  Eqn.  (34) 


with  Eqn.  (32)  and  Eqn.  (33)  we  have: 
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Hpm  mm 


ftpl  G1 
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Eqn.  (4o) 
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and  f .  =  f  / G 

tpg  tpin  I  _E 

°S. 


P 

m 


"tpl  / 

\G 


fk 


Eqn.  (41) 


g 


In  the  present  and  several  previous  investigations  at  the 
University  of  Alberta  the  following  approach  has  been  adopted.  The  well 
known  mechanical  energy  balance  (37)  may  be  written  for  each  phase  in 
the  following  way: 

WjL  dt  (v  dP  +  V^dV^  +  dF-^  =  0  Eqn.  (42) 


and 


t  « 


W  dt  (v  dP  +  V  dV  +  dF  )  =  0 

g  g  6 _ g  g 


Eqn.  (43) 


where  =  mass  rate  of  flow  of  the  liquid,  lbs.  per  sec. 


WCT  =  mass  rate  of  flow  of  the  gas,  lbs.  per  sec. 


g 


P 

v: 

v 


specific  volume  of  the  liquid  at  average  flow 
conditions,  cu.  ft.  per  lb. 

specific  volume  of  the  gas  at  average  flow  conditions, 
cu.  ft.  per  lb. 

pressure,  lb.  per  sq.  ft.,  absolute 

actual  linear  velocity  of  liquid,  ft.  per  sec.,  at 

average  flow  conditions 

actual  linear  velocity  of  gas,  ft.  per  sec.,  at  average 
flow  conditions 


dF^  =  irreversibilities  attributable  to  liquid  ft.  of  liquid 
dF0.  =  irreversibilities  attributable  to  gas  ft.  of  gas 
dt  =  differential  element  of  time,  sec. 
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Adding  Eqns.  (42)  and  (43)  yields: 


V1  dP  +  +  +  !fi  V  dP  +  ^  Vg  dVg  +  W£  aF  =  0  Eqn.  (44) 

SC  W1  W-,  ~ *7 


"I 

If  the  kinetic  energy  effects  are  neglected  or  considered  to  be  included 

within  the  terns  dF,  and  dF  we  have: 

1  g 


dP  (T1  +  Wg  vg)  +  (dF) 


=  o 


W- 


W 

where  (dF),  =  dF  +  g 

1  1  W-l  i 


Eon.  (45) 


where 


1  dP  (l  +  _g  g  )  +  (aF)  =  o 

u  ^r~ 

W1  'l 


-  VL  dP  (1  +  IV)  =  (dF). 


R 


W  v 

CT 


,  .  Q 


WI  vx 


g  ~  g  ,  air-water  volume  ratio. 


Ecn.  (46) 


If  the  right  hand  side  of  Eqn.  (46)  is  equated  to  the  Fanning 
equation  we  have: 

-  V]L  dP  (1  +  Rj  = 


2  f  V  dL 


gc3 


-  v 


dP  _ 


2  f  V 


1  dL  gcD  (1  +  By) 


Eqn.  (47) 


and  a  modified  "two-phase  friction  factor"  is  defined. 

Eqn.  (47)  is  integrated,  assuming  that  the  variations  in 
are  small  enough  to  justify  the  use  of  average  values  based  upon  an 
average  system  pressure  and  temperature.  Under  these  conditions: 
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AP 
V1  AL 


1  2  f  V2 

(1  +  Ry)  C  D 


Eqn.  (48) 


^tp  =  1 


AL  +  R  ) 


2fRVB  PR 
6c  D 


Eqn.  (49) 


The  subscript  'R'  refers  to  the  reference  phase  which  in  the 
present  investigation  is  the  liquid  phase.  Later  on,  subscript  ’s'  is 
used  to  signify  the  gas  or  second  phase  which  is  air. 

In  terras  of  the  mass  velocities  of  each  phase  Eqn.  (49)  can  be 
written  as: 

„  ->  o  -F 

R  R _  _ 

(1  +  Rv)  e 

C  S 


AP+.  1  2  f R  G  R„  2  f  G2  ,  . 

tP  =  7 - T  - =  7 - - - r  _ s.  ,  s  Eqn.  (50) 

AL  (l  +  R  )  g„  D  PR  (l  +  R^)  g  d  P 


v 


'c  R 

Comparison  of  Eqn.  (50)  with  Eqn.  (32)  and  Eqn.  (33)  shows  that: 


fj.  =  Rv  -p  = 

tpg  (iTTT  s 

v  V7 


/G 


(1  +  RV)\G 
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g  f- 
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Eqn.  (51) 


g  /  R 

A  review  of  the  various  equations  used  so  far,  shows  that  the 
following  relationship  exist  between  the  two -phase  friction  factor  f^g  , 
based  on  the  gas  phase  properties,  used  by  Schneider  (34)  and  those  used 


by  other  investigators. 
2 

f  =  ^  f  =  1 

tpg  T-  1 


(1  +  Rv)  \G 


"  i  '*  "”ft)  i  5; 


Schneider  Martinelli  Present  Investigation  Chenoweth  and  Kosterin,  Bertuzzi 

U.  of  A.  Martin  et  al 

Similar  relationships  can  be  written  for  the  friction  factor  based  on  the 


liquid  phase  properties: 

.2 


d  +  Rv) 


_  ~~  C  f 

R  ml 


in 


=  ftpm  lGnf  PA 

/  Pm 


-  36  - 


Hold-up  and  Slip  Velocity 

One  of  the  distinguishing  characteristics  of  two -phase  flow 
is  the  development,  within  the  test  section  of  a  concentration  of  one 
of  the  phases  greater  than  that  in  the  supply  mixture.  The  tens  ’’Hold-up 
ratio"  is  used  in  the  present  investigation  as  a  quantitative  measure 
of  the  liquid  accumulation.  Hold-up  ratio  is  defined  as  the  ratio  of 
a)  the  volume  ratio  in  the  supply  mixture  to  b)  the  volume  ratio  in 
the  flow  section.  Both  volume  ratios  are  expressed  in  cubic  feet  of 
gas  per  cubic  feet  of  liquid  at  average  flow  pressure  and  temperature. 
This  phenomenon  has  been  observed  by  many  investigators.  Martinelli 
(28,  29,  30)  has  measured  the  liquid  and  gas  ’hold-up'  and  has  found 
it  to  be  a  function  of  the  parameter  X  alone.  Alves  (2)  has  also 
treated  it  the  same  way.  Chisholm  and  Laird  (l2)  have  given  a  new 
terminology  to  the  same  phenomenon  and  have  called  it  ’saturation'. 

They  have  presented  an  equation  that  correlates  the  two -phase  pressure 
drop  and  'saturation'  and  have  developed  equations  for  the  prediction 
of  'liquid  saturation'.  Govier,  Radford  and  Dunn  (l9)  and  Govier  and 
Short  (20)  have  correlated  the  hold-up  ratio  as  a  function  of  the 
superficial  liquid  velocity  and  other  parameters.  Some  other  workers 
(19,  20,  2l)  have  referred  to  the  'slip  velocity'  of  the  gas  phase. 

Slip  velocity  is  defined  as  the  difference  in  the  actual  average  lineal 
velocity  of  the  gas  and  the  liquid  phases  both  measured  at  the  average 
flow  pressure  and  temperature.  The  two  quantities,  hold-up  ratio  and 
slip  velocity  are  two  measures  of  the  same  phenomenon  and  are  related 
as  follows. 


If  the  volume  fraction  of  the  flow  section  occupied  by  the 
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second  phase  (gas  phase)  he  then  (l-R  )  represents  the  volume  fraction 

s  s 

occupied  by  the  reference  phase  (liquid  phase). 


Hold-up  ratio  HR 


W  v 


S  g 

W6  V1 


l-R 

s 


1-r£ 

R, 


V 

s 


Eqn.  (52) 

Eqn.  (53) 


The  average  lineal  velocities  of  each  phase  are: 


v  =  Vs 
S  *1 


and 


VR  "  VR 


l-R. 


and  the  slip  velocity  V  is: 

si 


V  ,  =  V'  _  v 
si  s  VR 


R. 


V. 


R 

(T^) 


V- 


R 


l-R, 


V£ 

R~ 


(l-Rs)  -  l 


v. 
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V- 


R  ViiR 


l-R, 


(H„  "  1) 


Eqn.  (54) 


or  v;  -  =  v;  (hd  -  1) 


R  x  R 


Eqn.  (55) 
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EXPERIMENTAL  EQUIPMENT 

The  pilot  pipe  line  was  designed  to  obtain  pressure  drop  and 
hold-up  data  for  the  horizontal  flow  of  air  and  water  mixtures. 

A  schematic  flow  diagram  of  the  equipment  is  shown  in  Figure  1. 

The  main  test  section  of  transparent  cellulose  acetate  butyrate  pipe 
of  1.026  inch  inside  diameter,  with  an  overall  length  of  about  4l  feet 
was  supported  along  its  length  in  a  horizontal  position  by  angle  iron 
sections.  These  sections,  in  turn,  were  supported  by  steel  brackets  fixed 
onto  the  wall.  The  pipe  was  secured  in  position  by  clamps.  For  hold-up 
measurements,  two  round-port  quick  closing  plug  valves,  A  and  B,  connected 
by  a  mechanical  linkage,  were  installed  at  the  ends  of  the  plastic  pipe. 
Small  cocks  situated  in  the  valve  bodies  enabled  the  trapped  water  to 
be  drained  from  the  test  section.  The  plastic  pipe  was  attached  to  the 
round-port  valves  by  means  of  adapter  couplings.  The  discharge  end  of 
the  test  section  was  attached  to  a  separator  in  such  a  way  that  the 
flow  entered  the  separator  tangentially.  The  separator,  Plate  1,  was 
a  simple  hollow  cylindrical  vessel,  about  6  l/2  inch  in  diameter  and 
about  6  feet  high,  having  two  outlets,  one  at  the  bottom  for  water  and 
the  other  at  the  top  for  air.  The  outlet  for  water  was  provided  with  a 
valve  D  to  adjust  the  discharge  of  water.  A  "Belfield  parabolic  plug 
valve"  was  mounted  in  the  air  line  at  the  top  of  the  separator,  to  control 
the  mid  point  pressure  of  the  test  section. 

Two.,  2  inch  galvanized  iron  pipes  served  as  the  water  and  air 
lines.  Orifice  plates  were  installed  in  these  lines  to  measure  the  flow 
rate  of  the  two  components.  A  pressure  regulator  on  the  air  line  maintained 
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PLATE  I 


SEPARATOR 


the  pressure  of  the  incoming  air  between  50-55  psig.  Another  round-port 
plug  valve  C  was  installed  in  the  air  line  between  the  orifice  plate 
and  the  entrance  section,  to  minimize  the  fluctuations  in  the  manometers 
and  to  give  a  higher  pressure  upstream  of  the  air  orifice.  Govier, 
Radford  and  Dunn  (l9)  and  Hoogendoorn  (2l)  in  their  work  on  two-phase 
fluid  flow  concluded  that  the  method  of  introduction  of  the  two  phases 
did  not  effect  the  flow  pattern  provided  the  'calming  zone'  was  of 
sufficient  length.  Hence  a  simple  mixing  'tee'  was  used  in  this  work 
(Plate  2). 

A  sufficiently  long  'calming  zone'  was  provided  by  having  the 
upstream  pressure  tap  about  8  feet  from  the  upstream  valve  B.  The 
downstream  pressure  tap  was  about  2  feet  from  the  discharge  end  provid¬ 
ing  a  test  section  of  30  feet.  A  third  tap  was  made  halfway  along  the 
test  section  for  recording  the  mid-point  pressure.  The  pressure  taps 
were  1/8  inch  holes,  drilled  through  the  top  of  the  pipe.  Slipped  over 
these  taps  were  sleeves  sealed  with  'O'  rings. 

Air  and  Water  Supply 

The  air  was  supplied  from  the  University  of  Alberta  power 
plant  compressor  at  about  80  psig.  A  pressure  regulator  on  the  line 
maintained  the  pressure  of  the  air  between  50-55  psig. 

The  water  supply  was  obtained  from  a  constant  level  overhead 
tank  in  conjunction  with  a  5'H.P.  Fairbanks -Morse  centrifugal  pump, 
with  an  intake  pressure  of  25  psig.  and  a  discharge  pressure  of  about 
100  psig. 

Control  and  measurement  of  the  variables  were  obtained  with 
the  aid  of  orifice  plates  and  rotameters,  for  the  air  and  water  flow 
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PLATE  II  ENTRANCE  SECTION,  SHOWING  THE  MIXING 
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rates,  differential  water  and  mercury  manometers  for  the  pressure  drops 
and  simple  gravimetric  method  for  the  insitu  water  content  of  the  test 

section. 

Flow  Measurements 

Orifice  plates  were  used  to  measure  the  flow  rates  of  both  air 
and  water.  In  the  case  of  air,  one  rotameter  (Hoke.  Meter  No.  2102) 
and  six  orifices  ranging  in  diameter  from  0.082  inch  to  0.66  inch  and 
for  water  one  rotameter  (Fischer  and  Porter  S.No.  M6-I718/I)  and  five 
orifices  ranging  in  diameter  from  0.101  inch  to  0.752  inch  were  required 
to  cover  the  range  of  air  and  water  velocities  under  investigation. 

All  the  orifice  plates  were  fabricated  in  the  Department  Work  Shop. 
Flange  type  pressure  taps  were  used.  The  calibration  curves  for  the 
rotameters  and  the  orifice  plates  are  given  in  Appendix  C.  A  60  inch 
differential  water  manometer  and  a  60  inch  differential  mercury  mano¬ 
meter  recorded  the  pressure  drop  across  the  air  and  water  orifices 
respectively.  Air  was  the  transmitting  fluid  in  the  lead  lines  from 
the  air  orifice  plate  while  water  was  used  in  the  case  of  the  water 
orifice  plate. 

Pressure  Measurements 

The  pressure  drop  across  the  30  feet  test  section  was  regis¬ 
tered  on  either  a  60  inch  water  or  60  inch  mercury  manometer.  For  low 
pressures  an  inclined  water  manometer  was  used.  Since  air  was  the 
transmitting  fluid,  'splash  pots'  were  installed  at  each  of  the  pressure 
taps  to  prevent  any  water  from  being  carried  away  into  the  lead  lines. 
The  'splash  pots'  (Plate  3)  were  lucite  cylinders  about  2  inches  in 
diameter  and  6  inches  high.  A  needle  valve  connected  the  splash  pot 
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to  the  'O'  ring  sleeve.  Valves  were  provided  at  the  bottom  of  the 
splash  pots  to  drain  away  any  water  that  collected.  Lead  lines  of  l/4 
inch  transparent  polythene  tubing  connected  the  top  of  the  splash  pots 
to  the  manometers.  The  mid-point  pressure  was  registered  on  a  60  inch 
mercury  manometer  and  as  well  as  a  gauge.  A  'Brown  Recorder'  in 
connection  with  the  pressure  control  valve  maintained  the  mid-point  pres¬ 
sure  at  the  required  level.  Plate  4,  shows  the  arrangement  cn  the 
Rinel  Board  of  the  different  manometers,  the  two  rotameters,  the 
'Brown  Recorder'  and  the  valves  controlling  the  flow  rates  of  air  and  of 
water.  Figure  la  shows  the  manometer  connections. 

Temperature 

Two  dial  thermometers  (Bell  and  Gossett  30°F-250°F)  one  each  on 
the  air  and  water  lines  were  installed  to  record  the  temperature  of 
the  air  and  water.  A  third  one,  similar  to  those  above,  in  the  discharge 
end  recorded  the  temperature  of  the  mixture  as  it  left  the  test  section. 
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EXPERIMENTAL  PROCEDURE 


Before  conducting  the  two-phase  tests,  the  pressure  measuring 
equipment  and  the  rotameter  and  orifice  plates  calibrations  were  checked 
by  taking  pressure  drop  data  for  different  flow  rates  of  water  and  air 
alone.  During  those  test  runs  when  air  alone  was  used,  the  mid-point 
pressure  in  the  test  section  was  maintained  between  35-36  psia. 

Experiments  were  then  conducted  to  study  the  two -phase  flow  of 
air  and  water  over  a  range  of  superficial  water  velocities  from  0.01  ft. 
per  sec.  to  5*03  ft.  per  sec.  For  each  constant  superficial  water 
velocity  the  superficial  air  velocity  varied  from  0.1  ft.  per  sec.  to 
50  ft.  per  sec.  During  all  the  experiments  the  average  pressure,  or  the 
mid-point  pressure  in  the  test  section  was  maintained  at  35-36  psia.  This 
value  was  chosen  because  it  was  convenient  and  it  was  also  used  by  a 
earlier  worker  (38)*  For  each  test  pressure  drop  across  the  30  ft.  test 
section,  pressure  drop  across  the  air  and  water  orifices,  the  upstream 
pressure  at  the  air  orifice,  the  temperature  of  the  air,  water  and  the 
mixtures  and  the  insitu  water  content  were  recorded.  A  detailed 
description  of  one  such  experiment  follows. 

Orifices  to  give  the  desired  superficial  velocities  of  air 
and  water  were  first  selected  and  installed.  Valves  A,  B,  and  C  (Figure  l) 
were  set  in  the  open  position.  Valve  D  was  closed.  The  valves  on  the 
manometers  recording  the  pressure  drop  across  the  test  section  were  also 
closed.  Air  was  introduced  in  the  test  section  by  opening  the  valves 
P  and  Q.  Valve  C  was  then  adjusted  to  minimize  the  fluctuations 
particularly  in  the  manometer  across  the  air  orifice,  and  to  give  higher 
pressures  upstream  of  the  air  orifice.  When  the  mid-point  pressure 
reached  the  desired  level  (i.e.  35-36  psia),  water  was  introduced  in 
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the  system  by  turning  on  the  centrifugal  pump  and  opening  the  valves 
and  R.  0  D  was  tl.en  opened  and  adjusted  to  give  a  steady  discharge 

of  water  and  to  maintain  a  fairly  steady  level  in  the  separator. 

Valves  R  and  Q  were  adjusted  to  give  the  desired  flow  rates  of  reference 
phase  (water)  and  second  phase  (air).  When  the  conditions  inside  the 
test  section  had  reached  a  steady  state  as  indicated  by  visual  observations, 
a  note  of  the  flow  pattern  was  made.  The  time  interval  required  to 
establish  steady  state  conditions  varied  from  a  few  minutes  to  almost 

minutes.  In  general,  the  slower  the  flow  rates,  the  longer  the  time 
required  to  reach  steady  state.  In  addition  to  the  flow  pattern  the 
following  data  were  recorded: 

1)  Pressure  drop  across  the  water  and  air  orifices. 

2)  The  pressure  drop  across  the  30  ft.  test  section.  To  obtain  this,  the 
valves  on  the  manometers  were  opened  slightly.  The  pressure  differ¬ 
ential  as  registered  by  the  water  or  mercury  manometer  was  then 
recorded. 

3)  Mid-point  pressure. 

4)  Temperature  of  the  air,  water  and  the  mixture. 

5)  Upstream  static  pressure  at  the  air  orifice. 

6)  Hold-up:  After  recording  the  above  values,  and  checking  the  mano¬ 
meters  to  ascertain  that  the  correct  flow  rates  were  being  maintained, 
valves  A  and  B  were  closed  by  using  the  mechanical  linkage  and  the 
valves  Q  and  R  were  closed  immediately.  The  pump  was  stopped  and 

the  air  supply  turned  off  by  closing  the  valve  P.  The  small  cocks 
situated  in  the  body  of  valves  A  and  B  were  then  opened  slowly  and 
the  trapped  water  collected  in  a  bucket  and  weighed.  More  often,  the 
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pressure  inside  the  test  section  enabled  the  water  to  be  discharged 
easily.  When  necessary,  the  clamps  were  loosened  and  the  pipe  was 
slightly  lifted  at  the  center  to  facilitate  draining. 
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EXPERIMENTAL  RESULTS  AND  DISCUSSION 


The  results  of  the  preliminary  tests  of  single -phase  flow  of 
air  and  water  are  presented  in  the  form  of  the  conventional  friction 
factor  and  Reynolds  number  charts  in  Figures  5c  and  6c  of  Appendix  C. 

The  calculated  data  showed  good  agreement  with  theoretical  curves  in  the 


laminar  f  = 


16 

%e 


and  in  the  turbulent  region  with  the  Blasius  Equation, 


f  -  (0.079)  NRe  °’25  .  Because  of  the  good  agreement,  it  is  concluded 
that  the  calibrations  of  the  orifice  plates  and  the  pressure  measuring 
system  were  quite  accurate. 

The  flow  patterns  occurring  at  various  superficial  velocities 
of  second  phase  (air)  from  0.1  ft.  per  sec.  to  50  ft.  per  sec.  for 
constant  superficial  velocities  of  reference  phase  (water)  from  0.01  ft. 
per  sec.  to  5.03  ft.  per  sec.  were  noted.  Figures  2  to  11  indicate  the 
type  of  flow  pattern  observed  during  each  experiment.  In  these  figures 
the  letter  synbols  used,  represent  the  following  flow  patterns: 

B  =  Bubble 
PL  =  Plug 
ST  =  Stratified 
W  =  Wave  or  Ripple 
SL  =  Slug 

F  =  Film,  or  Semi annular  or  Annular 
Pressure  drops  and  liquid  hole-up  data  were  collected  for  the 
range  of  the  velocities  mentioned  above. 


Flow  Patterns 

Figures  2  to  11  indicate  the  type  of  flow  pattern  occurring  at 
various  superficial  velocities  of  second  phase  (air)  for  constant 
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superficial  velocities  of  reference  phase  (water)  from  0.01  ft.  per  sec. 
to  5.03  ft.  per  sec.  In  general  the  following  flow  patterns  were  observed 

l)  Bubble  2)  Plug  3)  Stratified  4)  Wave  or  Ripple 
5)  Slug  6)  Film  or  Semiannular  or  annular. 

The  sequence  of  appearance  of  the  flow  patterns  at  a  parti¬ 
cular  superficial  velocity  of  reference  phase  and  at  increasing  second 
phase  superficial  velocities  was  not  necessarily  in  the  order  given 
above.  It  was  not  possible  to  make  photographic  studies  of  the  flow 
pattern  during  each  experiment  but  some  representative  patterns  are 
given  in  Plates  5  and  6. 

Plate  5  shows  the  photographs  of  flow  patterns  developed  at 
superficial  velocity  of  reference  phase  of  0.3  ft.  per  sec.  and  super¬ 
ficial  velocities  of  second  phase  of  52.2  ft.  per  sec.,  21.0  ft.  per  sec 
and  1.2  ft.  per  sec.  The  flow  patterns  occurring  at  these  velocities  are 
’Film',  'Wave  or  Ripple'  and  'Stratified'  respectively.  Plate  6  shows 
the  flow  patterns  developed  at  superficial  velocity  of  reference  phase  of 
3.02  ft.  per  sec.  and  superficial  velocities  of  second  phase  of  10.0  ft. 
per  sec.,  2.00  ft.  per  sec  and  0.l4  ft.  per  sec.  The  flow  patterns  that 
developed  at  these  velocities  are  'Slug',  'plug'  and  'Bubble'  respectively 
The  'Slugs'  were  moving  at  such  high  velocities  that  it  was  not  possible 
to  get  a  good  photograph  with  the  equipment  available.  No  other  flow 
pattern  was  developed  at  second  phase  superficial  velocities  other  than 
those  mentioned  above. 

Bubble  Flow:  Bubble  flow  occurred  at  high  superficial  velocites 
of  reference  phase  and  low  superficial  velocities  of  second  phase.  This 
type  of  flow  occurred  at  superficial  velocities  of  reference  phase  of  more 
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Film  Flow 


=  52.2  ft.  per  sec. 


Wave  or  Ripple  Flow  Vg  =  21.0  ft.  per  sec. 


Stratified  Flow  V  =  1.2  ft.  per  sec. 

s 


PLATE  V 


FLOW  PATTERNS  AT  V  =  0.3  ft.  per  sec. 


\ 


Slug  Flow  V  =  10.0  ft.  per  sec. 

O 


Plug  Flow  V  =  2.00  ft.  per  sec. 

b 


Bubble  Flow  Vs  =  0.l4  ft.  per  sec. 


Bubble  Flow  =  0.l4  ft.  per  sec. 


PLATE  VI 


FLOW  PATTERNS  AT  L  =  3-02  ft.  per  sec. 
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than  1.0  ft.  per  sec.  and  second  phase  superficial  velocities  from  0.1  ft. 
per  sec.  to  1.0  ft.  per  sec. 

Plug  Flow:  This  type  of  flow  was  encountered  at  the  same 
superficial  velocities  of  reference  phase  as  in  the  case  of  bubble  flow 
but  at  higher  superficial  velocities  of  second  phase  i.e.  from  2.0  ft.  per 
sec  to  10.0  ft.  per  sec. 

Stratified  Flow:  Stratified  flow  generally  occurred  when  the 
superficial  velocities  of  the  two  phases  were  low.  For  the  range  of 
superficial  velocities  of  reference  phase  from  0.01  ft.  per  sec.  to  0.1  ft. 
per  sec.,  stratified  flow  occurred  when  the  second  phase  superficial 
velocity  range  was  from  0.1  ft.  per  sec.  to  5*0  ft.  per  sec.  For  super¬ 
ficial  velocities  of  reference  phase  between  0.1  ft.  per  sec.  and  1.0  ft. 
per  sec.,  the  maximum  superficial  velocity  of  second  phase  to  cause 
stratified  flow  was  2.0  ft.  per  sec.  At  higher  superficial  velocities 
of  reference  phase,  stratified  flow  was  not  encountered  during  the  range 
of  the  second  phase  superficial  velocities  employed.  Schneider  (3*0  has 
found  a  similar  behaviour. 

Wave  or  Ripple  Flow:  Wave  or  Ripple  flow  occurred  also  at  low 
superficial  velocities  of  reference  phase  i.e.  from  0.01  ft.  per  sec. to 
0.5  ft.  per  sec.  For  superficial  velocities  of  reference  phase  up  to  0.1 
ft.  per  sec.  wave  flow  occurred  when  the  superficial  velocity  of  second 
phase  was  between  5  ft.  per  sec.  and  10  ft.  per  sec.  But  for  superficial 
velocities  of  reference  phase  between  0.1  ft.  per  sec.  and  0.5  ft*  per  sec. 
superficial  velocities  of  second  phase  between  5  ft.  per  sec  and  20  ft.  per 
sec.  caused  wave  flow  pattern.  No  wave  flow  pattern  was  observed  for 
superficial  velocities  of  reference  phase  greater  than  0.5  ft.  per  sec. 
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at  any  of  the  superficial  velocities  of  second  phase. 

Slug  Flow :  This  type  of  flow  pattern  was  observed  only  at 
superficial  velocities  of  reference  phase  greater  than  0.5  ft.  per  sec. 
and  at  superficial  velocities  of  second  phase  greater  than  6  ft.  per 
sec.  Slug  flow  did  not  exist  for  any  superficial  velocities  of  reference 
phase  less  than  0.5  ft.  per  sec.  for  any  of  the  superficial  velocities 
of  second  phase. 

Film  Flow:  Film  or  Semiannular  or  annular  flow  occurred  at 
low  superficial  velocities  of  reference  phase  i.e.  between  0.01  ft.  per 
sec.  and  0.3  ft.  per  sec.  and  at  superficial  velocities  of  second  phase 
greater  than  10  ft.  per  sec.  Film  flow  was  not  encountered  for  superficial 
velocities  of  reference  phase  greater  than  0.3  ft.  per  sec.  even  at  the 
maximum  superficial  velocity  of  second  phase  employed  during  the  tests. 

A  chart  indicating  the  occurrence  of  the  various  flow  patterns 
is  given  in  Figure  12.  This  chart  has  been  prepared  on  the  basis  of  the 
experiments  and  observations  made  in  this  investigation  and  is  similar 
to  one  presented  by  Schneider  (3*0  for  the  natural  gas -kerosene  system. 

The  abscissa  and  the  ordinate  are  the  mass  velocities  (lbs.  per  hr.  - 
sq.  ft.)  of  water  and  air  respectively.  The  calculated  mass  velocities 
of  the  two  phases  and  the  corresponding  flow  pattern  are  presented  in 
Table  II  B  of  Appendix  B.  The  solid  lines  indicate  the  regions  developed 
on  the  basis  of  the  present  investigation.  The  dashed  lines  define  the 
same  regions  as  developed  by  Schneider  (3*0-  From  a  comparison  between 
the  two,  it  is  evident  that  the  shape  of  the  regions  is  similar  though 
not  identical.  The  flow  pattern  regions  developed  on  the  basis  of  the 
present  investigation  are  more  extensive  than  the  corresponding  regions 
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indicated  by  Schneider  (34).  'Cresting  Flow',  a  type  resembling  annular 
flow  but  occurring  at  much  higher  liquid  rates  and  mentioned  by  Schneider 
(3*0  in  his  plot  was  not  encountered  in  the  present  study.  At  the  same 
time  Bubble  Flow  is  not  mentioned  by  Schneider  (34)  as  his  studies  cover 
mainly  wave,  seraiannular  and  annular  flow. 

Pressure  Drop  and  Friction  Factor 

Two  phase  pressure  drops  were  recorded  for  various  superficial 
velocities  of  the  second  phase  (air)  at  constant  superficial  velocities 
of  the  reference  phase  (water)  of  0.01,  0.0303,  O.O697,  0.1,  0.3,  0.5, 

0.7,  1.0,  3*02  and  5*03  ft.  per  sec.  It  was  originally  planned  to  go  up 
to  a  reference  phase  superficial  velocity  of  10  ft.  per  sec.  but  it  was 
not  possible  since  beyond  a  superficial  velocity  of  5.0  ft.  per  sec.  the 
controller  could  not  maintain  the  mid-point  pressure  of  35 ~3&  psia. 

The  experimental  pressure  drops  are  presented  in  Table  II  A  of 
Appendix  A.  As  can  be  seen  from  this  table,  whenever  there  were  fluctu¬ 
ations  in  the  manometer,  which  could  not  be  avoided,  the  high,  the  low 
and  the  average  values  of  the  pressure  drops  are  given.  At  the  lower 
superficial  velocities  of  reference  phase,  low  second  phase  superficial 
velocities  could  not  be  employed  since  at  those  flow  rates,  the  pressure 
drops  were  too  small  to  be  measured  accurately  with  the  equipment  available. 

The  calculated  superficial  friction  factor  fr.  using  Eqn.  (49) 
is  tabulated  in  Table  I  B  of  Appendix  B,  along  with  other  calculated 
values  of  the  input  air-water  volume  ratio  R  and  the  hold-up  ratio,  E  . 

y 

The  superficial  friction  factor  f^  is  plotted  against  (R^  +  l)  in  Figures  2 


to  11,  each  at  a  different  constant  superficial  velocity  of  the  reference 
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phase.  The  addition  of  the  constant  1,  to  the  R  values  enabled  the  fD  vs 
(Rv  +  l)  curves  to  be  extrapolated  to  the  single  phase  friction  factor 
where  the  value  of  the  abscissa  is  one  since  this  value  of  abscissa 
represents  the  reference  phase  flowing  alone. 

Figures  2  to  5  represent  the  f  versus  (R^  +  l)  curves  for 
superficial  velocities  of  reference  phase  of  0.01,  0.0303,  0.0697  and 
0.1  ft.  per  sec.  At  these  low  superficial  velocities  of  reference  phase 
the  data  give  smooth  curves.  Figures  6  to  10  represent  the  same  plots  for 
superficial  velocities  of  reference  phase  of  0.3,  0.5,  0.7,  1.0  and  3*02 
ft.  per  sec.  Here  again  the  curves  are  smooth  but  exhibit  a  point  of 
inflection.  Figure  11,  representing  the  same  plot  for  a  superficial 
velocity  of  reference  phase  of  5*03  ft.  per  sec.  is  almost  a  straight 
line.  The  only  explanation  for  this  dissimilarity  is  that  at  low  super¬ 
ficial  velocities  of  reference  phase  the  flow  patterns  encountered  were 
'stratified',  'wave'  and  'film'  which  are  considered  to  be  more  stable  in 
nature  with  no  fluctuations  in  the  pressure  drops,  thereby  giving  smooth 
curves.  At  higher  superficial  velocities  of  the  reference  phase,  'plug* 
and  'slug'  type  of  flow  were  also  observed  though  there  were  no  instances 
of  'wave'  flow.  A  change  in  the  flow  pattern,  particularly  from  a  stable 
one  to  an  unstable  one,  might  have  caused  a  sudden  change  in  the  pressure 
drop.  Moreover,  a  look  at  the  plots  at  higher  superficial  velocities  of 
reference  phase  shows  that  a  change  in  the  shape  of  the  curve  coincided 
invariably  with  a  change  in  the  flow  pattern.  This  change  in  shape  with 
a  change  in  flow  pattern  is  more  pronounced  in  the  hold-up  ratio  curves 
which  are  discussed  later. 


From  data  shown  in  Figures  2  to  11,  values  of  the  superficial 
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friction  factor ,  f ,  were  obtained  for  input  air-water  volume  ratios 

of  1.0,  5*0,  10.0,  20.0,  50.0,  100  and  200.  Figure  13  shows  on  logarithmic 

co-ordinates,  the  variation  of  the  friction  factor,  f_,  with  the  Reynolds 

number  calculated  from  the  properties  of  the  reference  phase.  In 

addition  the  single  phase  friction  factor  relation  is  also  shown.  In 

effect,  Figure  13,  is  a  correlation  of  the  friction  factor,  f  ,  as  a 

R 

function  of  the  reference  phase  Reynolds  number  with  the  input  air-water 
volume  ratio  as  the  parameter.  The  data  of  Figure  13  in  the  laminar 
region  show  that  the  friction  factor,  f^,  is  inversely  proportional  to 
the  reference  phase  Reynolds  number  up  to  an  input  volume  ratio  of  50 
but  is  less  dependent  on  the  reference  phase  Reynolds  number  as  the  input 
volume  ratio  increases.  The  discontinuity  present  in  the  single  phase 
line  (Ry  =  0)  due  to  the  transition  from  laminar  to  turbulent  flow,  is 

more  pronounced  for  the  two -phase  lines  up  to  an  input  volume  ratio  of 
20.0.  As  the  input  volume  ratio  increases  the  discontinuity  decreases. 

In  the  turbulent  region  the  lines  for  the  input  volume  ratio  of  1.0  and 
5.0  are  very  similar  and  almost  parallel  to  the  single  phase  line.  At 
higher  values  of  input  volume  ratio  the  slopes  of  the  lines  decrease. 

For  an  input  volume  ratio  of  20.0  the  line  in  the  turbulent  region  is 
almost  parallel  to  the  abscissa.  For  higher  values  of  input  volume  ratio 
i.e.  50.0,  100  and  200  the  lines  have  a  positive  slope.  The  discontinuity, 
due  to  the  transition  from  laminar  to  turbulent  flow,  in  the  two-phase 
friction  factor  versus  Reynolds  number  curves  has  a  wider  range  than  in 
the  case  of  single  phase  flow.  Furthermore  the  transition  region  or  the 
discontinuity  appears  to  shift  to  the  left  gradually  with  increasing 
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values  of  input  volume  ratio.  No  explanation  is  forwarded  to  justify  the 
trend  noticeable  in  Figure  13  though  the  trend  seems  to  be  very  systematic . 

A  comparison  of  the  data  collected  in  the  present  investi¬ 
gation  with  the  Lockhart  and  Martinelli  (28)  correlation  is  presented 
in  Figures  l4  to  17 .  Calculated  values  of  the  parameters  X  and  Q> 
necessary  for  these  plots  are  tabulated  in  Table  III  of  Appendix  B. 

The  comparison  is  very  good  (Figure  l4)  when  both  the  phases  are  in 
turbulent  motion.  Figure  15  is  for  the  water  in  turbulent  and  air  in 
viscous  motion.  Here  also  the  data  compare  favourably  well  with  the 
Martinelli  correlation.  Figures  16  and  17  are  for  the  viscous -turbulent 
and  viscous-viscous  mechanism  respectively.  Here  the  comparison  is 
not  very  good,  almost  all  of  the  points  falling  lower  than  the  Lockhart 
and  Martinelli  curve.  However,  most  of  the  points  on  these  figures 
represent  experimental  data  in  the  stratified  region  of  flow  for  which 
Lockhart  and  Martinelli  correlation  is  not  valid  (l8). 

Comparison  of  the  superficial  friction  factor,  fR,  versus 
reference  phase  Reynolds  number  curves  for  the  present  investigation 
with  some  data  of  Schneider  (34),  Thomsen  (29)  and  Johnson  and  Abou- 
Sabe  (23)  is  presented  in  Figure  l8.  The  data  of  Schneider  (34), 

Thomsen  (29)  and  Johnson  and  Abou-Sabe  (23)  is  evaluated  on  the  basis 
of  the  calculations  employed  in  the  present  investigation.  No  data  is 
available  in  the  laminar  region.  Schneider's  (34)  data  though  in  the 
laminar  region,  corresponds  to  input  volume  ratio,  Ry,  values  greater 
than  500.  The  dashed  lines  plotted  in  Figure  l8,  indicate  interpolated 
data  at  lower  R  values  comparable  to  those  encountered  in  the  present 
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investigation.  The  influence  of  the  gas  phase  (natural  gas  in  Schneider's 
(3*0  work)  is  expected  to  he  negligible  and  that  of  the  liquid  phase 
(kerosene  in  both  Schneider  (3^)  and  Thomsen's  (29)  work)  is  taken 
into  account  by  the  abscissa  of  the  plot.  All  these  data,  fall  above 
that  for  the  present  investigation  probably  due  to  the  different  pipe 
diameters  employed.  The  work  of  Johnson  and  Abou-Sabe  (23)  is  the 
closest  to  the  present  investigation.  They  have  employed  the  same 
components  (water  and  air)  and  the  pipe  diameter  (O.O725  ft.)  is  also 
very  close  to  the  one  used  in  the  present  investigation.  No  explanation 
is  forwarded,  however,  for  the  different  trend  exhibited  by  their  data. 

On  the  whole  the  comparison  is  satisfactory. 

Hold-up  Ratio  and  Slip  Velocity 

As  mentioned  earlier,  hold-up  ratio  and  slip  velocity  are  two 
measures  of  the  same  phenomenon  which  results  in  the  development  within 
the  test  section  of  a  concentration  of  one  of  the  phases  greater  than 
that  in  the  supply  mixture.  Hold-up  ratio,  H^,  is  shown  as  a  function 
of  air-water  input  volume  ratio,  R  ,  in  Figures  2  to  11,  for  different 
constant  superficial  velocities  of  reference  phase.  Up  to  a  reference 
phase  superficial  velocity  of  0.3  ft.  per  sec.  the  curves  are  smooth 
as  shown  in  Figures  2  to  6.  At  reference  phase  superficial  velocities 
greater  than  0.3  ft.  per  sec.  however  these  curves  exhibit  a  discontinuity. 
This  change  which  is  shown  in  Figures  7  to  11  is  more  pronounced  in  these 
curves  than  in  those  for  the  superficial  friction  factor.  At  lower 
reference  phase  superficial  velocities  only  '  stratified',  'wave',  and 
'film'  flow  were  encountered  while  at  higher  superficial  velocities  of 
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the  reference  phase,  'plug'  and  'slug'  flow  also  were  observed.  Here 
again,  a  change  in  the  shape  of  the  curve  has  always  coincided  with  a 
change  in  the  flow  pattern.  This  tends  to  show  that  the  hold-up 
ratio  or  liquid  hold-up  is  also  a  function  of  the  flow  pattern.  A 
similar  trend  is  noted  by  Levy  (27).  There  is  a  strong  tendency  for 
all  these  curves  to  extend  to  a  hold-up  ratio  of  unity  at  zero  input 
volume  ratio.  This  can  be  expected  since  in  the  limiting  condition 
of  zero  Ry  or  all  liquid  flow,  the  slip  velocity  is  zero  and  by  Eqn.  (55) 
the  hold-up  ratio  equals  unity.  The  same  trend  was  noted  by  Govier, 
Radford  and  Dunn  (19)  in  their  studies  on  the  vertical  upward  flow  of 
air-water  mixtures. 

A  composite  plot  of  hold-up  ratio,  H^,  versus  R^  at  different 

reference  phase  superficial  velocities  was  first  prepared  but  there 

was  so  much  overlapping  of  the  curves,  that  it  is  not  included  here. 

Figure  19  is  a  composite  plot  of  H  versus  the  second  phase  superficial 

R 

velocity  V  with  the  superficial  velocity  of  reference  phase,  V  ,  as 
S  R 

the  parameter.  The  curves  are  separate  and  distinct,  thereby  emphasizing 
the  importance  of  the  parameter.  Once  again  in  the  limiting  condition 
when  Vo  is  zero  or  all  liquid  flow/,  the  curves  approach  a  hold-up  ratio 
of  unity. 

Slip  velocity,  ,  is  calculated  by  using  Eqn.  (5*0-  The 

si 

calculated  values  are  presented  in  Table  IV  B  of  Appendix  B.  The  slip 
velocity,  V  ,  is  plotted  in  Figure  20  against  air-water  input  volume 

S  «L 

ratio,  R  ,  with  the  superficial  velocity  of  reference  phase,  VR,  as  the 
parameter.  At  higher  values  of  R^_  most  of  these  curves  have  a  slope  of 
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about  unity.  A  general  trend  towards  zero  slip  velocity  at  an  input 
air -water  volume  ratio  of  zero  is  also  noticeable  from  Figure  20. 

This  is  quite  logical  since  at  a  of  zero,  that  is  all  liquid  flow, 
the  slip  velocity  is  also  zero. 
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CONCLUSIONS 

1)  Six  types  of  flow  pattern  namely,  bubble,  plug, 
stratified,  wave  or  ripple,  slug  and  film  (semiannular  or  annular)  were 
noted.  A  tentative  flow  pattern  chart  is  presented  for  air -water 
system.  For  such  a  combination  of  fluids,  it  allows  the  prediction 

of  the  type  of  flow  which  would  exist  under  a  given  set  of  mass 
velocities  of  the  two  phases. 

2)  Pressure  drop  data  for  the  horizontal  flow  of  air-water 
mixtures  have  been  obtained  and  correlated  by  plotting  a  superficial 
friction  factor,  based  on  the  water  properties  against  a  superficial 
reference  phase  (water)  Reynolds  number.  This  correlation  is  limited 

to  air -water  systems,  in  a  one  inch  smooth  horizontal  pipe  at  an 
average  system  pressure  of  35“36  psia  and  temperature  of  55°-75°F.  The 
data  compares  quite  favourably  with  the  Lockhart  and  Martinelli 
correlation  (28),  although  deviating  somewhat  in  the  viscous-turbulent 
and  viscous-viscous  mechanisms. 

3)  Hold-up  ratio  is  plotted  against  second  phase  super¬ 
ficial  velocity  and  the  input  volume  ratio.  In  both  cases,  in  the 
limiting  conditions  of  the  second  phase  superficial  velocity  or  the 
input  volume  ratio  being  zero,  the  hold-up  ratio  approaches  unity. 

k)  Slip  velcoity  is  plotted  against  input  volume  ratio. 

In  the  limiting  case  of  all  liquid  flow,  the  slip  velocity  approaches 
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NOMENCLATURE* 

A  -  Cross-sectional  area  of  the  pipe, 

C  -  Volume  fraction  of  a  component  in  the  supply 
mixture  at  average  flow  conditions, 

D  -  Inside  diameter  of  pipe, 

f  -  Friction  factor  in  general, 

gc  -  Dimensional  conversion  constant, 


G  -  Mass  velocity  of  a  component, 

Hp  -  Hold-up  ratio,  ratio  of  input  volume  ratio  to  insitu 
volume  ratio,  both  volumes  at  average  flow 
conditions, 

AL  -  Length  of  pipe  line, 

NRe  -  Reynolds  number, 

AP  -  Pressure  drop, 

AP 

XT  -  Pressure  drop  per  unit  length, 

Q  -  Volume  rate  of  flow  at  average  flow  conditions, 

R  -  Fraction  of  pipe  volume  occupied  by  a  component 
at  average  flow  conditions, 

Rj4  -  Input  gas-liquid  mass  ratio, 

Rv  -  Input  gas-liquid  volume  ratio,  both  volumes  at 
average  flow  conditions, 

v  -  Specific  volume  of  the  fluid  at  average  flow 
conditions, 

V  -  Superficial  velocity  of  fluid  at  average  flow 
conditions  based  on  total  pipe  cross-section, 

V'  -  Actual  average  lineal  velocity  of  fluid  at 
average  flow  conditions, 

Vsi-  Actual  average  slip  velocity  =  V('r  -  Vj  , 

W  -  Mass  rate  of  flow, 

X  -  Lockhart  and  Martinelli  parameter, 


sq.  ft. 


dimensionless 

ft. 

dimensionless 
lb.  (mass)  ft. 
lb.  (force)  sec/" 

lbs. per  sec.  sq.ft. 

dimensionless 

ft. 

dimensionless 
lbs. per  sq.  ft. 
lbs.  per  sq.ft.  -  ft 
cu.  ft.  per  sec. 

dimensionless 

dimensionless 

dimensionless 


cu. 

ft. 

per 

ft. 

per 

sec. 

ft. 

per 

sec. 

ft. 

per 

sec. 

lbs 

.  per 

sec 

dimensionless 


* 


Special  terms  have  been  defined  in  their  proner  context. 
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Lockhart  and  Martinelli  parameter, 


p  -  Density  of  fluid  at  average  flow  conditions, 


p  -  Absolute  viscosity  of  fluid  at  average  flow 
conditions, 

y  -  Kinematic  viscosity  of  fluid  at  average  flow 
conditions, 


dimensionless 
lbs.  per  cu.  ft. 

lbs. per  ft.  sec. 

sq.  ft.  per  sec. 


Subscripts: 

g  -  Refers  to  the  gas  phase  or  values  calculated  by  using  the  gas  phase 
properties. 

1  -  Refers  to  the  liquid  phase  or  values  calculated  by  using  the  liquid 

phase  properties. 

m  -  Refers  to  the  two-phase  mixture  or  values  calculated  by  using  the 
properties  of  the  mixture. 

R  -  Denotes  the  reference  phase  (Water,  in  the  present  investigation) 

or  values  calculated  by  using  the  properties  of  the  reference  phase. 

s  -  Denotes  the  second  phase  (Air,  in  the  present  investigation)  or 
values  calculated  by  using  the  properties  of  the  second  phase. 

tp  -  Refers  to  two -phase  in  general. 

tpg  -  Refers  to  two-phase  but  using  the  properties  of  the  gas. 

tpl  -  Refers  to  two-phase  but  using  the  properties  of  the  liquid. 

tpm  -  Refers  to  two -phase  and  using  the  properties  of  the  mixture. 
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Table  I  A  gives  general  information  about  the  pipeline 
test  section. 

Table  II  A  gives  the  experimental  data  collected  for  the 
flow  of  air-water  mixtures  in  a  horizontal  pipe. 

For  a  particular  experiment  number  the  table  gives: 

The  superficial  velocities  of  water  and  air,  ft.  per  sec. 

the  temperature  of  water,  air  and  the  mixture,  °F. 

the  pressure  drop  over  the  test  section  of  30  feet,  lbs.  per  sq.ft. 

the  mid -point  pressure,  lbs.  per  sq.  in.  abs. 

the  flov;  pattern,  and 

the  weight  of  water  trapped,  in  the  test  section,  lbs. 
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TABLE  I  A 
GENERAL  DATA 

Overall  length  of  pipe  line,  4l.O  ft. 

Test  section  length,  i.e.  distance  over  which  the  pressure 
drop  was  measured,  30*0  ft. 

Inside  diameter  of  pipe,  1.026  inch  =  0.0855  ft. 

Volume  of  pipe  line  for  hold-up  measurements,  0.2345  cu.  ft. 
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Appendix  B  gives  all  the  calculated  data  under  different 
tables.  Each  table  is  followed  by  a  sample  calculation  as  an 
illustration. 

Table  I  B  gives,  besides  the  superficial  velocities  of  water 
and  air,  the  following  calculated  data: 

Input  volume  ratio,  (R  ).  . 

v 7  input 

Insitu  volume  ratio,  (Rv)insitu 

Hold-up  ratio, 

Density  of  the  second  phase, 

i.e.  air  at  average  flow  condition,  p_,  lbs.  per  cu.  ft. 

The  superficial  friction  factor,  fR 
(These  data  are  plotted  in  Figures  2  to  11  and  19 . ) 

Table  II  B  gives  the  mass  velocities  of  the  two  phases  and 
the  flow  pattern.  (This  data  is  employed  in  Figure  12. ) 

Table  III  B  primarily  gives  the  calculated  values  of  the 
parameters  *..4\  and  X.  It  also  tabulates  the  Reynolds  numbers  of 

O 

the  two  phases,  the  pressure  drop  in  lbs.  per  sq.  ft.  for  the  two  phases 
flowing  alone  and  the  observed  two-phase  pressure  drop  in  lbs.  per  sq.  ft. 
(The  se  data  are  used  in  presenting  the  comparisons  in  Figures  l4  to  17 . ) 

Table  IV  B  gives  the  values  of  the  slip  velocity  and  the  input 
volume  ratio  which  are  used  in  Figure  20. 
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SAMPLE  CALCULATION 


Experiment  No.  112 


Input  volume  ratio  R 
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Insitu  volume  ratio  =  0.227 6  =  33*00 

0.00&9 

(p  ) .  ^ 
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ZIP  =  76*90  lbs.  per  sa .  ft.  from  Table  IX  A 

tp 


&L  =  30.0  ft. 

D  =  0.0855  ft. 
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•  R 
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SAMPLE  CALCULATION 
Experiment  No.  112 

G  =  p  V  =  (62.4)  (0.0303)  (3600) 

R  R  R 

=  6800  lbs.  per  hr.  sq.  ft. 

G  x  10  3  =6.80  lbs.  per  hr.  sq.  ft. 

Gs  =  °S  vs  -  (0.181*)  (51.6)  (3600) 

=  3l*j200  lbs.  per  hr.  sq.  ft. 

G^  x  10  =  34.2  lbs.  per  hr.  sq.  ft. 
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SAMPLE  CALCULATION 

Reynolds  number  of  the  reference  phase 

 D  VR  ^ 


(Vr 


R 


(0.0835)  (0-0303)  (62.4)  =  2h0 

(I)  (0.000672) 


Pressure  drop  for  waiter  flowing  alone 


(APR) 


2  f  AL  V  PR 

gc  D 


f  at  (Op  Of  240  =  16 


l6  =  0.066 


240 


AP  =  (2)  (0.066)  (30)  (0.0303)^  (62.4) 

(32.17)  (0.0855) 

=  0.08175  lbs.  per  sq.  ft. 

Reynolds  number  of  the  second  phase 

(N  )  =  D  V  p 

Re  s  s  s 


=  (0.0855)  (51.6)  (0.184)  =  69,000 

(0.0175)  (O.OOO672) 

At  (lv4>e)s  °f  69^000  friction  factor  from  Blasius  eqn. 
f  =  0.079  .  =  0.004874 

W70S0) °*25 

AP  =  (2)  (0.004874)  (30)  (51. 6)2  (0.184)  =  52.05  lbs. 


(32.17)  (0.0855) 

1/2 


ft, 


t  - 

AP  /AL 
tp' 

1/2 

76.9/3O 

'  g 

aps/al  __ 

52.05/30, 

=  1.21 


<b 

AVAL' 

1/2 

76.9/30  1 

1 

apr/al 

0.08175/30 

1/2 


=  30.70 
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TABLE  IV  B 


SLIP  VELOCITY 


1 

2 

3 

4 

5 

Expt 

.  Superficial* 

Hold-up 

* 

Volume 

* 

Slip 

No. 

velocity  of 

ratio 

fraction 

velocity 

reference 

phase 

V  ft. /sec. 

R 

“r 

of  pipe 
occupied 
by  second 
phase  Ro 

V  ft./; 
si 

101 

0.01 

39.4 

O.726 

1.40 

102 

0.01 

40.0 

0.731 

1.43 

103 

0.01 

64.2 

0.764 

2.68 

104 

0.01 

73.5 

O.873 

5.71 

105 

0.01 

59.6 

0.944 

10.64 

106 

0.0303 

13.62 

0.556 

0.862 

107 

0.0303 

23.75 

O.586 

1.662 

108 

0.0303 

32.80 

O.676 

2-975 

109 

0.0303 

38.70 

0.8127 

6.10 

110 

0.0303 

32.60 

0.908 

10.40 

111 

0.0303 

38.3 

O.9473 

21.42 

112 

0.0303 

51.5 

O.9706 

52.00 

113 

O.O697 

14.55 

0.522 

1-975 

114 

O.O697 

23.50 

0.572 

3.50 

115 

O.O697 

24.20 

O.7675 

6.95 

116 

O.O697 

22.80 

0.884 

13.10 

117 

O.O697 

22.60 

0.943 

26.4 

118 

O.O697 

31.00 

0.961 

53.6 

119 

0.1 

7.28 

0.415 

1.073 

120 

0.1 

13.30 

0.437 

2.19 

121 

0.1 

17.70 

0.5415 

3.64 

122 

0.1 

19.65 

O.7292 

6.89 

123 

0.1 

17.50 

0.855 

11.40 

124 

0.1 

19.68 

0.9153 

22.02 

125 

0.1 

25.90 

0.9515 

51.4 

126 

0.3 

2.175 

O.2515 

0.461 

127 
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0.3 

3.39 

0.336 

1.081 

0.3 

6.95 

0.349 

2.742 

129 

0.3 

7.25 

0.4945 

3.715 
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0-3 

6.08 

0.7315 

5.67 

131 

0.3 

9-7^ 

0.7875 

12.32 

132 

0.3 

11.08 

0.8627 

22.02 

133 

0.3 

I6.3O 

0.9173 

55.5 

*  At 

average  flow  conditions 

6 

Input 

volume 

ratio 
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209.0 
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991.0 

17.1 

33-6 

68.6 

167.6 

328.0 

686.0 

1,700.0 

14.9 

29-55 

74.03 

149.80 

298.50 
724 . 0 

5.19 

10.28 

20.95 

53-00 

103.00 

212.50 
506.0 

0.727 

1.715 

3.72 

7.06 

16.55 

36.00 

69.50 
181.00 


V  /  .  . 


•  I 
» 


1  ■  - 


. ,  .. 

)  v  - 


I 


0: 


I  I 


I 


I 


.  —  l  .... 

L  >".LO\ 

-  xo 

-  ■•-■• 

fO-. -J  .  -L. 

-  v>  .  J 

Xv-'-  • 

.  o. 

cl  j 

. 

v  *  ^ 

i  * 

lv.. 

--  i  • 

U  *  J  . 

*  o 

V  _ 

* 

_  m  V/ 

C;  " 

•  -  l  "  - 

1  "  ^ 

'n.  -  1 

•01 

*  U_*_ 

' 

• 

Vv-J. 

1  C  .  C 

•  • 

C.L 

•  -- 

1  * 

•  ■ 

I  ^  -A- 

^  1 

r  • 

>  , 

• 

O.  ■ 

^  a  O 

i  * 

.  - 

'  -  - 

c .  ..... 

. 

-  .  J 

vJJ _ 

•  - 

1  1  .  * 

•  {/ 

- - 

C  .  .  .  ^ 

1 

.... 

1  )  -  - 

'  ^  * 

0.  .  * 

i 

vj  •  J 

C 

I  \  * 

0  .  . 

i’-. 

— 

G  1  1  ’ 

• 

i  ^ 

S.  '  * 

...  . 

Gv  - 

. 

w.  - 

•  0 

*  -  •  _  j 

.  - 

• 

l  -  — 

o  G 

l  ... 

1  . 

•  . 

.x.X 

i'  • 

v  *  • 

. 

r 

o .  • 

V-  ■  ^  • 

X.  .  . 

i  • 

^  .  X 

.  •  ' 

■  .... 

v  • 

*•  .  •  V- 

X  .  C 

• 

•  v. 

.  0 

. *  C- 

•  . 

C  'v  . 

CY  -  • 

• 

* 

V 

1  • 

. 

^Q. . 

.  * 

C  • 

C  •  V 

. 

1  *  'n 

C  i  * 

• 

• 

^  1  1  • 1 

1  • 

V  •' 

• 

V*  •  CC 

\  ■  • 

• 

. 

-  115  - 


TABLE  IV  B  (Cont'd.) 


1 

2 

3 

4 

5 

6 

134 

0.5 

1.20 

0.181 

0.1221 

0.29 

135 

0.5 

1-585 

0.212 

0.371 

0.424 

136 

0.5 

1.918 

0.3545 

0.711 

1.044 

137 

0.5 

2.68 

0.421 

1.45 

1.952 

138 

0.5 

2.68 

0.6092 

2.15 

4.184 

139 

0.5 

3.98 

0.7125 

5.18 

9.460 

140 

0.5 

- 

- 

- 

20.24 

l4l 

0.5 

6.72 

0.7592 

11.88 

21.20 

142 

0.5 

8.775 

0 . 8277 

22.60 

42.15 

143 

0.5 

12.08 

O.8962 

53-40 

104.00 

144 

0.7 

2.24 

0.39 

1.422 

1-43 

145 

0.7 

2.32 

O.562 

2.082 

2.98 

146 

0.7 

3.965 

0.652 

5.96 

7-^3 

147 

0.7 

5.60 

O.726 

11.75 

14.82 

148 

0.7 

7.310 

0.801 

22.2 

29.40 

149 

0.7 

10.00 

O.879 

52.05 

72.60 

150 

1.0 

1.42 

0.115 

0.1605 

0.1468 

151 

1.0 

1.158 

0.155 

O.I87 

0.210 

152 

1.0 

1.452 

0.265 

0.615 

0.520 

153 

1.0 

1.85 

0.3645 

1.338 

1.055 

154 

1.0 

1.635 

0.5282 

1.345 

1.825 

155 

1.0 

2.860 

0.6480 

5.285 

5.25 

156 

1.0 

4.615 

0.6850 

11.46 

10.00 

157 

1.0 

6.13 

0.7680 

22.10 

20.22 

158 

1.0 

8.7O 

0.8560 

53-50 

51.50 

159 

3.02 

1.037 

0.0448 

0.117 

0.0486 

160 

3.02 

1.774 

0.044 

2.442 

O.079 

161 

3.02 

1.02 

0.145 

0.0707 

O.173 

162 

3.02 

1.435 

O.1925 

1.628 

0.343 

163 

3.02 

1.467 

0-335 

2.122 

0.734 

164 

3.02 

1.62 

O.508 

3.805 

1.668 

165 

3.02 

2.04 

0.6l4 

8.15 

3.24 

1 66 

3.02 

4.235 

0.7935 

47.25 

16.25 

167 

5.03 

1.438 

0.067 

2.36 

0.1018 

168 

5.03 

1.432 

0.128 

2.49 

0.210 

169 

5-03 

1.332 

0.24 

2.195 

0.430 

170 

5.03 

1.49 

0.3985 

4.10 

0.982 

171 

5.03 

1-73 

0.530 

7.81 

1.945 
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SAMPLE  CALCULATIONS 
Experiment  No .  112 

From  calculations  following  Table  I  B,  we  have 

Hold-up  ratio  =  51-5  and 

volume  of  pipe  occupied  by  gas  =  O.2276  cu.  ft. 

volume  fraction  of  pipe  occupied  by  gas  Rs  =  O.2276  = 

0.2345 


From  Eqn .  ( 5^ ) 

Slip  velocity,  Vgl 


Slip  velocity  V£l 


VR  (Hr  -  1) 

(1  -  fss> 

(0.0303)  (51-5  -  1) 

(1  -  O.97O6) 

52.00  ft.  per  sec. 
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Calibration  of  Rotameters  and  Orifices 

The  rotameter  and  orifice  plates  for  air  up  to  0.125  inch 
in  diameter  were  calibrated  by  using  a  gas-holder.  The  orifice  plates 
for  air  from  0.188  inch  to  0.66  inch  in  diameter  -were  calibrated  by 
using  a  critical  flov  prover.  The  original  calibration  data  are 
tabulated  in  the  research  note -book.  The  calibration  curve  for  the 

rotameter  is  given  in  Figure  lc,  in  which  the  air  flow  Q  in 

cu.  ft.  per  min.  (at  60°F  and  l4.4  psia)  and  the  flow  rate  divided  by 
the  rotameter  reading  is  plotted  against  the  rotameter  reading  as  the 
abscissa.  The  calibration  curves  for  air  orifice  plates  is  given  in 

Figure  2c,  in  which  the  flow  rate  Q  in  cu.  ft.  per  min.  (at  60°F 

O 

and  l4.4  psia)  is  plotted  against  PAh  on  logarithmic  co-ordinates 

T 

a 

where , 

P  is  the  static  pressure,  upstream  of  the  orifice,  psia 
Ah  is  the  pressure  drop  across  the  orifice,  inches  of  water 
Ta  temperature  of  flowing  air,  °R 

The  rotameter  and  the  orifice  plates  for  water  were  calibrated 
by  weight-time  measurements.  The  original  data  is  tabulated  in  the 
research  note -book.  The  calibration  curve  for  the  rotameter  appears  in 
Figure  3c  in  which  the  flow  rate  in  cu.  ft.  per  sec.  divided  by  the 
rotameter  reading  is  plotted  against  the  rotameter  reading.  The  cali¬ 
bration  curves  for  the  orifices  appear  in  Figure  4c  as  a  plot  of  rate 
of  flow  of  water  in  cu.ft.  per  sec.  versus  the  pressure  drop  across  the 
orifice  in  inches  of  mercury  on  logarithmic  co-ordinates. 
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Pipe  Diameter 

The  internal  diameter  of  the  pipe  was  determined  by  inside 
calipers.  A  large  number  of  readings  were  averaged  to  give  an  inside 
pipe  diameter  of  1.026  inches. 

Pipe  Volume  for  Hold-up  Measurements 

The  total  volume  between  the  two  plug  valves  A  and  B  was 
found  to  be  0.234-5  cu.  ft.  A  number  of  single  phase  runs  with  water 
were  conducted  and  the  water  trapped  was  collected  and  weighed.  The 
difference  between  the  theoretical  and  the  average  of  a  number  of  amounts 
of  water  collected  was  0.21  lbs.  or  1.435$  of  the  theoretical  value.  The 
loss  was  attributed  to  the  water  sticking  to  the  pipe  walls  and  the 
valves.  This  value  (0.21  lbs.)  was  added  to  the  weight  of  water  collected 
in  each  of  the  two  phase  flow  experiments  to  measure  the  hold-up. 

Pipe-line  Calibration 

The  calibration  of  the  flow  meters,  the  pipe -line  diameter 
and  the  operation  of  the  pressure  drop  measuring  equipment  was  checked 
by  conducting  single  phase  runs  with  air  and  water  flowing  alone.  The 
pressure  drop  and  the  calculated  friction  factors  are  tabulated  in  the 
research  note-book.  They  are  presented  graphically  in  Figures  5c  and  6c 
together  with  the  curve  f  =  l6  for  the  laminar  region  of  flow  and  the 
Blasius  equation  for  the  turbulent  region.  Good  agreement  is  evident 
between  the  experimental  data  and  accepted  curves. 
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where : 


P  «  Upstream  static  pressure,  p.s.i.a. 

Ah  *  Pressure  drop  across  the  orifice,  in  H20 
T0  *  Temperature  (*R)  of  flowing  air 
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